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Gade [Nucl. Phys. B 398, 499 (1993)] has shown that the density of states for a particle hopping on a two- 
dimensional bipartite lattice in the presence of weak disorder and in the absence of time-reversal symmetry 
(chiral unitary universality class) is anomalous in the vicinity of the band center e = whenever the disorder 
preserves the sublattice symmetry. More precisely, using a nonlinear a model that encodes the sublattice (chi- 
ral) symmetry and the absence of time-reversal symmetry she argues that the disorder average density of states 
diverges as | e | 1 exp (— e| lne| K ) with c some nonuniversal positive constant and K = 1/2 a universal exponent. 
Her analysis has been extended to the case when time-reversal symmetry is present (chiral orthogonal univer- 
sality class) for which the same exponent K = 1/2 was predicted. Motrunich et al. [Phys. Rev. B 65, 064206 
(2002)] have argued that the exponent K = 1/2 does not apply to the density of states in the chiral orthogonal 
universality class. They predict that K = 2/3 instead. We confirm the analysis of Motrunich et al. within a 
field theory for two flavors of Dirac fermions subjected to two types of weak uncorrelated random potentials: a 
purely imaginary vector potential and a complex valued mass potential. This model is the naive continuum limit 
of a model describing a particle hopping on a square lattice in the background of a Jt-flux phase and subjected 
to weak disorder that preserves the sublattice symmetry and time-reversal invariance. By commonly held uni- 
versality arguments, this model is believed to belong to the chiral orthogonal universality class. Our calculation 
relies in an essential way on the existence of infinitely many local composite operators with negative anomalous 
scaling dimensions. 

PACS numbers: 71.30.+h, 72.15.Rn, 64.60.Fr, 05.40.-a 



I. INTRODUCTION 

The effect of disorder on the density of states (DOS) in the 
problem of Anderson localization is known to be rather mild: 
In general, disorder smoothes the DOSji^ For example, the 
logarithmic Van Hove singularity in the DOS at the band cen- 
ter for a single-particle undergoing uniform nearest-neighbor 
hopping on a square lattice is washed out by disorder in two 
dimensions (2D)»M Consequently, the signature of disorder- 
induced critical behavior in Anderson localization is usually 
not captured by the DOS, and one needs to resort to more 
complex probes such as conductance statistics, energy levels 
statistics, or multifractal spectra that are sensitive to a dereal- 
ization transition. 

Exceptions to the rule that the DOS is smooth in Anderson 
localization can occur at the band edges or at the band center^ 
The band center is very special whenever the disorder pre- 
serves "particle-hole" symmetry, i.e., a symmetry that guaran- 
tees that energy eigenvalues always occur in pairs of opposite 
signs for any given realization of the disorder. The particle- 
hole symmetry can take two different forms: (i) the sublattice 
symmetry^i^iSiSiiS (also called chiral symmetry) which is rel- 
evant in the context of Hamiltonians defined on bipartite lat- 
tices and is implemented by a unitary transformation (chang- 
ing the sign of the wave functions on every site belonging to 



one sublattice but not on the other) and (ii) the Bogoliubov-de 
Gennes (BdG) symmetryiiiiSii^ which is relevant in the con- 
text of disordered superconductors and is implemented by an 
antiunitary transformation.— The special role played by the 
band center when the disorder preserves particle-hole symme- 
try follows from the enhanced level repulsion felt by energy 
eigenvalues close to the band center as a result of their mirror 
images. Naively, one would expect that level repulsion results 
in a depletion of the DOS close to the band center. This expec- 
tation is in fact always realized in finite systems and is relevant 
to quantum dots or to level statistics in the core of vortices in a 
superconductor^ However, it has been known since the work 
of Dy sonA on disordered one-dimensional chains that the DOS 
can, in spite of the enhanced level repulsion, diverge at the 
band center in the thermodynamic limit. The band center must 
then represent "some sort" of disorder-induced critical point 
characterized by a derealization transition of "some kind." 
Indeed, imagine that there is no derealization transition, i.e., 
that the localization length is finite at the band center. If so, 
the naive expectation of a depletion of the DOS as a result 
of enhanced level repulsion must apply in the thermodynamic 
limit as the system behaves in effect as a sampling of uncorre- 
lated boxes of linear size of the order of the localization length 
for sufficiently long length scales. 

Anderson localization with particle-hole symmetry thus 
represents a laboratory in which disorder-induced critical be- 
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havior can be studied both in ID and 2D. In ID a wealth 
of complementary approaches based on the nonlinear a 
model (NLSM) ) 16 i 17 real-space renormalization group (RG) 
techniques^ extensions of the Dorokhov-Mello-Pereyra- 
Kumar (DMPK) equationfI2J(i2£ direct studies of critical zero 
modes t 2 -' 22 and others^^iS^ are available and have been used 
to provide a rather detailed description of the derealization 
transition at the band center. In 2D, the use of the NLSM 
was pioneered by Gade and Wegner in the context of Ander- 
son localization with particle-hole symmetry implemented by 
a sublattice symmetry. § They predicted that the disorder aver- 
age DOS per unit volume and energy diverges upon approach- 
ing the band center e = as 

|e|^exp(-c|lne| K ), D = l, k=-. (1.1) 

Here, c is a positive nonuniversal constant, whereas the ex- 
ponents v and K are universal. Moreover, neither v nor 
K depends on whether time-reversal symmetry is broken or 
not. The NLSM approach of Gade and Wegner has been ex- 
tended to the random 7t-flux phas o 26 i 27 and to the random flux 
modelSlt also has been refined to account for the expectation 
value of the staggered DOS (Ref. l29i) and for the presence of 
Wess-Zumino-Witten (WZW) termsS In all cases, the DOS 
II. 1> was recovered. Particle-hole symmetry is also essential 
to understand the role played by disorder in a mean-field and 
non-self-consistent treatment of superconductivity^ In 2D, 
the disorder average DOS was predictediii^ to diverge as 

|lne| K (1.2) 

when spin-rotation symmetry is broken. The universal expo- 
nent K = 1/2 in the presence of time-reversal symmetry and 
K = 1 in the absence of time-reversal symmetry. 

Motrunich et al. in Ref. |30| have proposed that the exponent 
K is not given by 1 /2 as suggested by various fie ld-theoretica l 
calculations leading to Eq. ( fTTTl (see Refs.lil and f26l27l28l29l) 
but by the larger exponent 

k 4 (L3) 

The method by which Motrunich et al. reach their conclu- 
sion in Ref.|30jis based on (i) a real-space RG analysis that is 
intrinsically a strong-disorder technique^ (ii) a combination 
of exact variational bounds for strong random vector poten- 
tial disorder with the computation of the DOS in the weak- 
disorder limit by Ludwig et al. in Ref. H and with the compu- 
tation of the p function for the variance of the random vector 
potential by Guruswamy et al. in Ref. l27l and (iii) a connec- 
tion to the problem of monomer free energies in some prob- 
lems of dimer covering. In this paper we undertake the task 
of understanding the discrepancy between Eqs. J 1 . 1 1 and II. 3> 
for the Hatsugai-Wen-Kohmoto (HWK) model^ in the limit 
of weak disorder at microscopic length scales and within a 
purely field-theoretical approach. We shall show that the RG 
analysis for the average DOS in the HWK model performed 
in Ref. that leads to Eq. II. Il l is not complete in that it fails 
to capture the broadness of the local density of states (LDOS) 



distribution. A consistent RG analysis for the DOS can, how- 
ever, be constructed and yields Eq. dl .31 in the HWK model. 

The HWK model is the naive continuum limit for a model 
of a particle undergoing nearest-neighbor hopping on a rect- 
angular lattice in the background of 7t-flux phase. The mi- 
croscopic disorder is introduced by allowing real-valued ran- 
dom fluctuations of the nearest-neighbor hopping amplitudes 
about the 7t-flux phase. The continuum limit of the clean spec- 
trum is approximated by two flavors of Dirac fermions. The 
microscopic disorder is assumed to be sufficiently weak for 
the nodal structure of the clean energy spectrum to survive 
the perturbation by the disorder. It was shown in Ref. |32| 
that the disorder manifests itself by a purely imaginary ran- 
dom vector potential and by a random complex valued mass 
in the naive continuum limit. Furthermore, boundary condi- 
tions were chosen so as to guarantee the existence of a zero 
mode for any realization of the disorder and a numerical study 
of zero modes confirmed the theoretical expectation that zero 
modes are critical (multifractal). Numerical studies of the 
DOS for the HWK model are inconclusive due to severe lim- 
its on the achieved energy resolutions^! From a purely sym- 
metry point of view the HWK model belongs to the chiral 
orthogonal universality class since time-reversal symmetry is 
preserved. If symmetries alone control the asymptotic behav- 
ior of the DOS in the vicinity of the band center, one would 
expect that the disorder average DOS is of the form given in 
Eq. II. U . This is indeed what is found in Ref. |^ where it 
is shown that (i) the HWK model realizes a nearly confor- 
mal field theory and (ii) the DOS dl.lt follows from solving 
a Callan-Symanzik equation that depends on three coupling 
constants. 

In this paper we shall build on the remarkable results (i) of 
Guruswamy et al. in Ref. |^ to identify an infinite but count- 
able set of relevant local operators with negative anomalous 
scaling dimensions in the HWK model. The physical inter- 
pretation of these relevant operators is that they govern the 
dependence on energy e of all moments of the LDOS pro- 
vided e is sufficiently close to the band center e = 0. Their 
relevance reflects the fact that the probability distribution for 
the LDOS becomes very broadly distributed in the thermo- 
dynamic limit as its distribution is governed by the (multi- 
fractal) critical state at the band center. Correspondingly, 
we shall show that the disorder average DOS v(e) obeys a 
Callan-Symanzik equation that depends on infinitely many 
relevant coupling constants. Following an approximation de- 
veloped by Carpentier and Le Doussal for the random phase 
XY model in 2Dp4 we reduce the Callan-Symanzik equation 
with infinitely many relevant coupling constants to the Callan- 
Symanzik equation of Guruswamy et al. with the caveat that 
the anomalous scaling dimension in (ii) that corresponds to 
the annealed average value of the dynamical exponent must 
be replaced by the anomalous scaling dimension that corre- 
sponds to the quenched average value of the dynamical expo- 
nent. In doing so we recover Eq. II. 3> . 

The physical origin of the infinite but countable family of 
operators with negative anomalous scaling dimensions can be 
traced to the multifractality of zero modes supported by the 
HWK model when the random complex mass is switched off. 
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The explicit interplay between multifractality of zero modes 
and negative scaling dimensions in a conformal fie ld theory 
(CFT) was studied in Refs. I9I and 13 6l3 7l3 8l3 9l40l4Tl for a sin- 
gle flavor of a Dirac fermion subjected to a random vector 
potential. It was shown in Refs. 38 39 40) that the multifractal 
spectrum is characterized by a freezing transition. The mul- 
tifractal spectrum probes the scaling of the support of a nor- 
malized zero mode whose squared amplitude (height) is some 
prescribed value as a function of L in a finite box of volume 
L 2 . As the height is increased, the multifractal spectrum un- 
dergoes a freezing transition. Below the freezing transition, 
the support scales with the system size to some finite power. 
Above the freezing transition, the support does not scale any- 
more with the system size. A similar freezing transition oc- 
curs when the height is decreased towards a minimum value 
of th e sq uared amplitude of zero modes. It was argued in 
Refs. l30l andl42l that the freezing transition in the multifractal 
spectrum mirrors itself in the DOS for a single Dirac fermion 
subjected to a random vector potential. The relation between 
the multifractal spectrum and the scaling of the DOS is em- 
bodied by the dynamical exponent that relates scaling with 
energy to scaling with system size. It is thus the dynamical 
exponent that undergoes a freezing transition. We will argue 
that the same mechanism applies, with suitable modifications, 
to the HWK model. 

Although the emphasis is on Anderson localization, 
the models considered in this paper are also intimately 
connected to disordered Coulomb gases, random XY, 
and other 2D disordered systems from classical statisti- 
cal physicsi2iy2i2£i22idii22i22i2i Perhaps the most surprising 
connection between Anderson localization and a classical 
stochastic model is the one conjectured in Ref. |38| between 
the multifractal spectrum of zero modes and the freezing tran- 
sition of directed polymers on a Cayley tree (see Refs. l45l 
and l46h. Carpentier and Le Doussal have convincingly ar- 
gued that a consistent RG analysis of the multifractal spectrum 
of zero modes reduces to solving the so-called Kolmogorov- 
Petrovsky-Piscounov (KPP) nonlinear diffusion equation. 44 
Their approach yields new insights into the universality of the 
multifractal spectrum in terms of the unique properties of the 
asymptotic solutions of the KPP equation. It also gives a RG 
justification for the conjecture in Ref. l38l bv building on the 
fact that it is known from Ref. |45| that the transformation law 
of the generating function for the random partition function of 
directed polymers on a Cayley tree under an infinitesimal RG 
step is also governed by the KPP equation. Hereto we shall 
see that the Callan-Symanzik equation obeyed by the DOS in 
the HWK model is intimately related to the KPP equation. 

Our paper is organized as follows. Lattice and contin- 
uum versions of the HWK Hamiltonians are defined in Sec. 
im where we also define the Motrunich-Damle-Huse (MDH) 
model. We choose a supersymmetric (SUSY) field theory to 
represent the DOS in Sec. [Hi] After these preliminaries, we 
derive the DOS for the MDH model in Sec.|IV] We compute 
the DOS for the HWK model in Sec.|V] We conclude in Sec. 
IVII Some technical details are summarized in the Appendixes. 



II. DEFINITIONS 
A. Lattice models 

We are concerned with the naive continuum limit of a sub- 
class of lattice Hamiltonians belonging to the chiral orthogo- 
nal universality class. Consider some two-dimensional bipar- 
tite lattice — say, a lattice that is rectangular in the bulk. One 
sublattice is denoted A, the other B. There are Na sites on 
sublattice A and Nb sites on sublattice B. Without loss of gen- 
erality, assume that Na > Nb- The random Hamiltonian is then 
represented by 

where Tab is a Na X Nb real-valued matrix with the random 
matrix elements (Tab)!] = Uj if U are nearest-neighbor sites 
with ; e A and (Tab),) = otherwise, while the Nb X Na matrix 
Tba is the transpose of Tab- Tba = (Tab) 1 - Random Hamilto- 
nians of the form J2. It are generic representatives of the chi- 
ral orthogonal universality class. Wave functions in the N- 
dimensional Hilbert space (N := Na +Nb) are represented by 

By construction, H changes by a sign under the unitary trans- 
formation 

Va^+Va, ^b^-^b- (2.3) 
Zero modes have the form 

¥ = ( ^ J , (2.4) 

where 

G = Tba x ¥a = {Tab? x ¥a. (2.5) 

Zero modes are thus eigenstates of the unitary transformation 
d2.3i . There are Na — Nb zero modes. 

The HWK lattice model is given by the Hamiltonian #~ H wk 
of the form ( 12. H whereby it is assumed that (i) a 7t-flux phase 
threads every elementary plaquette P of the rectangular lattice 
and (ii) there are small random fluctuations 8f, 7 - of the nearest- 
neighbor hopping amplitude: 

nf=-< 4 - (2-6) 

ijeP 

Here, t is the characteristic energy scale (°< bandwidth) for 
the clean system. In Ref. the multifractal spectrum of 
zero modes of #hwk was studied numerically and shown to 
agree within numerical e rror bars wit h a theoretical predic- 
tion made in Refs. and l37l38l39l40l In subsequent papers 
the DOS was also studied numericallyi^i Unfortunately, the 
lack of energy resolution was so severe that the observation of 
a diverging DOS could only be established in Ref. 34 through 
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a tiny upturn of the DOS away from an algebraic decaying 
crossover regime. 

Motrunich et al. in Ref. [3(J proposed a different 
parametrization of the random fluctuations of the nearest- 
neighbor hopping amplitudes. Their lattice Hamiltonian ^ MDH 
is of the form (12. li but with 



f.. = e +M Vfj 
1.1 ij ' 



(2.7) 



ijeP 



instead of Eq. (12. 6> . The (j)'s are real-valued, white-noise- 
correlated, Gaussian random variables of vanishing means. 
With the help of the Na x Na diagonal matrix 



U A :-- 



(2.8) 



and of the Nb X A^b diagonal matrix 



V 



(2.9) 



the MDH random Hamiltonian becomes 




B 



o 



-1 



(2.10) 



where Tj^j is a Na x A^b matrix made of f^ '. In this model, it 
is straightforward to construct explicitly zero modes. To each 
solution 



of 



AB y T . 



(2.11) 



(2.12) 



[i.e., Eq. d2.5i without disorder], there corresponds a zero 
mode of given by 



(^) 



-l 





U n 







(2.13) 



This explicit construction of zero modes will have a coun- 
terpart in the field theory later on and plays a very special 
role. The imaginary gauge chiral transformation ( 12. 1 3I > turns 
the Schrodinger equation at any finite e into 



T<$% = e(U A r 2 V A , 
(T^)\' A = e(U B ) +2 % 



(2.14) 



B. Naive continuum limit 

Common to the HWK and MDH lattice Hamiltonians in the 
clean limit is the property that the Fermi surface collapses to 
four Fermi points at (±7t/2,±7t/2) in the Brillouin zone at 
half filling. Correspondingly, the clean DOS has a V-shaped 
singularity at the band center. Linearization of the spectrum 
around those four Fermi points describes two copies of "rel- 
ativistic" particles moving with the Fermi velocity vy := 2at, 
a being the lattice spacing. These are the two flavors of Dirac 
fermions. Unless specified, we will set v/ to 1. 

The difference in the parametrization of the disorder in the 
HWK and MDH lattice Hamiltonians survives the continuum 
limit. For sufficiently weak disorder, the naive continuum 
limit for the HWK lattice Hamiltonian reads^ 



with 



and 



Dl 







/OoA + a 3 A3 



(2.15a) 



(2.15b) 



(2.15c) 



Pauli matrices cr = (01,02), together with the unit 2x2 ma- 
trix Oo and the third Pauli matrices 03 have been introduced. 
Summation over the repeated indices (n = 1,2) is assumed 
throughout this paper. This Hamiltonian is invariant under 
time reversal, i.e., 



T(//„ WK )*T =//„ 



T 



01 

-Oi 



It also changes by a sign under charge conjugation 

Oo 



C "HWK C — "HWK 3 C 



-00 



(2.16) 



(2.17) 



which is nothing but the continuum counterpart of the trans- 
formation that implements the transformation Tab — ► — Tab in 
Eq. 12.31) . The o's encode the spinor grading of a single fla- 
vor of Dirac fermions. The (real-valued) disorder is repre- 
sented by (i) the "purely imaginary" random vector potential 
A = {A\,A'i) e M 2 , (ii) the "purely imaginary" random scalar 
potential (chemical potential) Aq 6 R, and (iii) the random 
mass potential A3 G R. 

The terminology "purely imaginary" comes from the fac- 
tor of i multiplying A and Ao in the operator £> HWK ; i.e., A 
and Ao parametrize the anti-Hermitian contributions to £> HWK . 
The Hermitian contributions to £> IIWK are the kinetic energy 
icr ■ d and the random mass contributions O3A3. With a 
prejudice towards a field-theoretical approach, it is standard 
practice to assume that A, Ao, and A3 are all independent 
white-noise-correlated Gaussian random variables of vanish- 
ing means4i The implicit assumption made here is that the 
detailed shape of the disorder should not matter as long as it 



is weak and short-range correlated in space. So far the dis- 
cussion is generic to the chiral orthogonal universality class. 
We now specialize to a subclass of the chiral orthogonal uni- 
versality class by, following Hatsugai et ai, making a further 
important simplification. Motivated by analytical continua- 
tion of a one-loop RG calculation performed in Ref. for 
a single flavor of Dirac fermions subjected to "purely real" 
random scalar, vector, and mass potentials^ we assume that, 
on the one hand, the two components of the vector potential 
have the same variance g A and, on the other hand, the random 
scalar potential and the random mass have the same variance 
g M (disorder averaging is denoted by an overbar): 

A 8^8{x-y), fi,V = 1,2, 



A jU (x)A v (y) = i g M 8^8(x-y), ^ = 0,3, (2.18) 



0, 



otherwise. 



This completes the definition of the HWK model. 

For sufficiently weak disorder, the naive continuum limit 
for the MDH Hamiltonian is given by Eq. 12 . 1 5 al without the 
purely imaginary random scalar potential and random mass. 
To see this, we demand a counterpart to Eqs. (12 . 1 Oi and (12 . 1 3 i 
in the continuum limit. First, we introduce the parametrization 

A M = d p <Pi +d M <t> 2 (9^ := e^v, /u = l,2) (2.19) 

of the random vector potential in terms of two independent 
scalar fields <I>i and <I>2, where £u V is the antisymmetric tensor 
with £12 = 1. Second, we observe that 

iO M (d M +A M ) = e +to3*i-ao*2 ia ^e- ia ^ + ^, 
ia^ —Aft) = e -*>3*i+°o<J>2 /a^e +iO3<1> i- O0<I> 2. 



We thus identify Eqs. (|28} and (|29} with 



(2.20) 



(2.21) 



To sum up, in the naive continuum limit the MDH model is 
given by 



D 



D 

t 



MLMi 





, D ltam = ia ll (d ll +A ll ), (2.22) 



where the vector potential A is randomly distributed according 
to 



A fl (x)A v (y)=g A 8 fJV 8(x-y), /i,v = l,2. (2.23) 



III. DOS AND FIELD THEORY 

A. Generalities 

The DOS v(e) per unit energy and per unit volume L 2 for 
the single-particle Hamiltonian H with eigenvalues e„ is de- 
fined by 



1 



v(e ) := 72L 8 ( e - e ») 
1 



L- 
1 



— V lim — Im 

L 2 Vl^ 07t 



e - e„ - ;r| 



lim — =■ y - 

n^onl 2 V (£-£„) 



(3.1) 



With T| kept small but finite, Eq. ( 13. Il l defines the smeared 
DOS. The LDOS v(e;x) is defined in terms of the exact nor- 
malized eigenfunctions of H by 



v{e;x) := lim -> --^ — \ 

n^onV (e-e«) 2 +il 2 



(3.2) 



When the single-particle Hamiltonian H is one member of a 
statistical ensemble, it is commonly held that the smeared 
DOS is self-averaging in the thermodynamic limit. On the 
other hand, critical statistical fluctuations of the wave func- 
tion amplitudes at x open up the possibility that the (smeared) 
LDOS is broadly distributed in the thermodynamic limit. 

The computation of the (smeared) LDOS in field theory is 
very natural as it reduces to the calculation of the expectation 
value of a local composite operator. For example, the smeared 
LDOS can be equivalently well represented by fermionic or 
bosonic path integrals 



1 f ©[\(/,\|/]i|/(a;)i|/(a;) exp \- fd 2 r\\i(r) (H-e + ir\)\\i(r)] 

— Im '■ — '■ 

Jt /2>[\pr,\|i] exp [- / d 2 r\\f(r) (H — e + ir\)i\r(r)] 

1 ^ J g> [p, P] V(xMx) exp [- Jd 2 r p(r) (g-e + h\) p(r)] 

Jt / ©[p,p] exp[-/^ 2 rp(r)(i/-£ + /ri)p(r)] ' 



respectively. Here, \\t and \\t denote two independent Grassmann fields whereas p denotes the complex conjugate of the complex- 
valued field p. 
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A SUSY representation for the smeared DOS can be achieved by combining the first with the second line of Eq. (13 .31 : 



v(e-/r|) = -^Im|©[¥,V,P,P] J £ / 2 a;¥ ( a; )^( a; ) e -/^'-[v('-)(«-e+^)v('-)+P(-)(«-^)P(-)] 



~2^2 Im ( / d2x [v( x M x ) + P(*)P(aO] 



(3.4a) 



Zsusy(£-'T|) 



where ((• • • ))z SVSY (£-ir\) denotes averaging with the SUSY partition functions defined by 

Z SUSY (e-jr|) := / ©[¥,¥,P,P] expi - f d 2 r [vj/(r) (H -E + ir\)y(r) + P(r) (H -E + ft\)$(r)] 



(3.4b) 



r 



Alternatively, one can replicate either line of Eq. i3.3\ to ob- 
tain a fermion (boson) replica representation of the DOS. 

By assuming that disorder averaging restores translation in- 
variance, we infer that the smeared LDOS v(e — iry,x) and 
the smeared DOS v(e — ;'r|) share the same mean value as sta- 
tistical variables. However, their distributions are likely to 
be very different in the thermodynamic limit L — ► °° at some 
fixed energy e — ;T| as the former is expected to be broadly 
distributed whereas the latter is expected to be self-averaging. 
From now on, we will refer to the disorder average DOS as 
the DOS having in mind the expectation that the DOS is self- 
averaging. 



B. SUSY partition function for the LDOS in the HWK model 

Needed is a SUSY partition function from which the LDOS 
in the HWK model can be represented by the expectation 
value of some local composite operator. To construct a con- 
venient SUSY partition function we choose a realization of 
the disorder Ao, A\, A%, and A3 in Eq. ( 12. 151 and define the 
fermionic partition function 



7(f) .. 



if) 



2)[¥4,¥4,¥b,¥s 



- fd 2 r\L if) +L lf) +L 
]g I [ AB BA +i 



{/)] 



L AB - VA^hwkVB: 



4a =¥b£>Lk¥a, 



(/) 



-e(Wa + W«). 



(3.5) 



In view of the particle-hole symmetry, we will assume that 
e > from now on and without loss of generality. The SUSY 
partition function is obtained by first, multiplying Z^' by Z^ h \ 
whereby Z^' follows from Z^O after replacing all Grassmann 
fields ¥ and \\t by the complex-valued fields p and p, re- 
spectively, and second, by integrating over disorder. We are 
now going to perform three successive changes of path inte- 
gral variables in the "matter" (xj/'s, \|/'s, p's, and p's) sector. 
We will only specify the first two changes of variables in the 
fermion (¥> sector as it is understood that the same changes 
of variables apply to the boson (P) sector. We will treat the 



fermion and boson sectors in an asymmetric way in the last 
step. Although these are symbolic manipulations, they can be 
justified by choosing appropriate unitary equivalent represen- 
tations of H UWK for the fermion and boson partition functions, 
respectivelyiS 
First, introduce 



y A =:v(/ A Oi, 



Vb--Vb<*u 



(3.6) 



whereby \\>a and \|/g are left unchanged. Denote eigenstates of 
G3 according to 

owa± = T¥a±, OiVb± = T¥b±, (3.7) 
say, and rewrite the disorder as 

A±:=Ai±/A 2 , M±:=A ±;A 3 . (3.8) 
This gives the fermion Lagrangians 



r(f) 

J AB 



rif) 

-BA 



(3 + +A + )v B _+v* + (d_+A_)y B+ ] 
+i(y*A-M + y B+ +wX + M_yif B _) , 

i i d + - A +) ¥a-+¥b+ (3- -A_) W A +] 
- i (¥b- M_m? a+ + Yb+ m + ¥a- ) . 



-e[¥1-¥a+ + ¥b-¥b+ + (- 



-)], 



(3.9) 



whereby d± : = d \ ± id2- Their boson counterparts follow from 
substitution of the \|/*'s by the P*'s and of the ¥ s by the p's. 
Second, define 



¥a : 
¥a 



1 



¥i, ¥b 



1 



¥2, 



2ni V2k 

1 1 
¥2, ¥b =: -=¥1 



271/ 



2n 



(3.10) 



The rescaling of the fields is done for later convenience when 
computing operator product expansion (OPE). We also rescale 
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the energy e — > 2m. This gives the fermion Lagrangians 



Af) 

^AB ~ 



Af) 

^BA ~ 



^ [Vi- (3+ ¥i-+¥t+ (3- +A_ 



2jc 



(¥l-^+¥i++¥t+^-¥i 



— (3+-A+) ¥2- +¥2+ (3- -A_) ¥ 2+ ] 

- ^ (¥2- M_ ¥2+ + ¥2+ M + ¥2-) . 



(/) 



e (¥l-¥ 2 + - ¥l-¥i + + ¥t+¥ 2 - - ¥2+¥i 



(3.11) 



Their boson counterparts follow from substitution of the V|/* 's 
by the p*'s and of the ¥ s by the P' s - Observe that to go from 



.(/),, (P) 



■ (f) 



X^'g to ig^' + L B b jl it suffices to change by a sign the 



J AB 



vector potential A and the scalar potential Aq, while leaving 
the sign of the mass potential A3 unchanged. 

This asymmetry between the two sublattices can be re- 
moved by performing the final change of path integral vari- 
ables according to2i 



¥2 

PI 



+¥ 2 > 
+P2, 



¥2 
P2- 



+¥2> 

-PI. 



(3.12) 



with all other fields unchanged. This gives the fermion La- 
grangian 



Af) 

J AB 



-(/) 
-BA 



E (3++A + )¥ fl - 



a=\2 



-Ya+ (3_+A_)v (H . 
-\? a _M + y a+ +\f a+ M_M? a _], 

(3.13a) 



The "decoupling transformation" ( 12.201 now reads (the fla- 
vor index runs over a = 1,2) 



in the fermion sector, and 



¥ fl - 
¥ QH 



hr'<J>i — <I>2 



P«- 



(3.14a) 



(3.14b) 



in the boson sector. Under the chiral transformation (13. 14l > we 
note that the fermion Lagrangians become 

z;l a=l,2 v 

M + <r 2i< V«+ + ¥ ;;+^^ +2,<I>i ¥:,- 

4 /5 = W 2 ** (¥'i- ¥2+ + Vr+¥2-) 

-ee- 2 *i(V 2 -Vi++V2+Vi-), (3A5a) 

whereas the boson Lagrangians become 



a=l,2 



4 6) =- £e + 2 *2 (PtP2* + + Pf + P2-) 
- £e - 2<I> 2(P' 2 _P' 1+ + P' 2+ P^ 



(3.15b) 



together with its boson counterparts obtained by substitution 
of the \|/*'s by the p*'s and of the ¥ s by the p's: 

^+42 = E ^-[PUa + +A + )p fl _ 

+ p;, + (a_+A_)p a+ 
+p:_M + p a+ + p; + M_p u _]. 

(3.13b) 

The asymmetry between sublattices is not completely re- 
moved, however, since it is disguised as a perturbation for the 
DOS with non-Hermitian appearencei^ 

4 f) = e(¥t-¥2+-¥2-¥i+ + ¥t+¥2--¥2+¥i-), 

4 h) = -e(Pt_PI + + p 2 _p 1+ + p 1+ PI_ + p 2+ p 1 _). 

(3.13c) 



The complex-valued random masses M± couple to the 
transversal component <t>i of the random vector potential 
through a bounded phase factor exp(±2/«J>i ) whereas they de- 
couple from the longitudinal component <t>2 of the random 
vector potential. Conversely, the energy term couples to the 
longitudinal component <t>2 of the random vector potential 
through an unbounded exponential factor exp(±2<t>2) whereas 
it decouples from the transversal component <t>i of the random 
vector potential. We conclude, with the help of the decou- 
pling transformation ( I3.14K that the random vector potential 

enters in the energy perturbation Xg + £g as it does in the 
squared amplitude of the zero mode of the operator £>„ DH in 
Eq. 12.221 — namely, through the unbounded exponential fac- 
tors exp(±2«J>2). We shall see that these exponential factors 
dictate the dependence on e of the DOS in the close vicinity 
of the band center^ 

The SUSY partition function for the HWK model that we 
will choose to work with treats quantum mechanical (through 
the SUSY matter fields \|/ and p) and disorder (through the ran- 
dom vector and mass potentials) averaging on an equal foot- 
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ing: 



IV. DOS FOR H M 



-Jd 2 rL 



2gA 
+L 



-A+A- 



2g 



-M+M- 



M 



if) 



■ (f) , L (f) 



(b) 



■(b) 



lb) 



where L 



(f/b)^ L (f/b)^ r (f/b) 



AB 



BA 



i E are given by Eq. 



(3.16) 



Corre- 



spondingly, the SUSY partition function for the MDH model 
is obtained from Eq. ( I3.16> by switching off the random 
masses M±, i.e., setting the variance g M = 0. This SUSY par- 
tition function is critical at the band center e = as follows 
from the decoupling transformation 13.141 . Indeed, it is then 
given by 



L, = ad 2 a + r^-(3„*i) 2 + ^(d„*2) 2 
z 8a z 8a 



Af), r (/) .,(*)., 00 

' AB ^BA ' AB ' ^BA > 



(3.17) 



whe re L$ + L$ and L^j + L%> follow from Eqs. J3.15al 
and J3. 15bi with the random masses M± switched off, respec- 
tively. We have exponentiated the Jacobian that arises from 
expressing the vector potential in terms of its longitudinal and 
transversal components with the help of fermionic ghosts a 
and a. In doing so we ensure that the SUSY partition function 
remains unity. The Fujikawa Jacobian arising from the decou- 
pling transformation in the fermionic sector cancels with that 
from the bosonic sector. 

It was shown in Ref. |^| that Eq. J3. 17i defines a CFT for 
each value of the marginal coupling constant g A . The full op- 
erator content was constructed explicitly; i.e., anomalous scal- 
ing dimensions of all local operators that can be built out of the 
\|/'s and P's were calculated. Furthermore, it was established 
that this CFT realizes a gl(2|2) current algebra. This prop- 
erty allows the explicit construction of primary fields realiz- 
ing one- and two-dimensional representations, say, of the cur- 
rent algebra. The anomalous scaling dimensions of primary 
fields carrying multiples of the U(l) gauge charge of 2 were 
shown to be unbounded from below, related to the "Cooper- 
like" terms entering L^', and responsible for the multifractal- 
ity of the zero mode. We are going to combine these results 
with a method pioneered by Carpentier and Le Doussal in Ref. 
l35l to compute the DOS for the MDH model. We will then 
show that the anomalous scaling dimensions of the family of 
Cooper-like terms that govern both multifractality of the zero 
mode and the DOS in the MDH model also control the DOS 
in the HWK model. To this end, we will make use of the fact, 
demonstrated by Guruswamy et al. in Ref.|^, that Eq. (13 . 1 6i 
defines a nearly conformal field theory in which the coupling 
constant g A is marginally relevant (in the sense that it grows 
logarithmically with length) and the coupling constant g M is 
exactly marginal. 



(/) 



(b) , r (b) 



In this section, we are going to compute the DOS in the 
MDH model. We shall develop an approximation that, pre- 
sumably, becomes exact in the limit e — > and in which the 
DOS for the MDH model reduces to the DOS for a single 
species of Dirac fermion subjected to a real-valued random 
vector potential. The cornerstone of our calculation is to dis- 
tinguish between the annealed and quenched average values 
of the dynamical exponent in very much the same way as one 
must distinguish between the scaling with system size of the 
average and typical values of the normalization of zero modes. 
In this way we will have established a one-to-one correspon- 
dence between the power-law dependence of the DOS and the 
multifractal spectrum of zero modes. 

Our starting point is the critical theory d3 . 1 71 from which 
we extract the nonvanishing two-point functions 



(*2(Z,2)*2(0)>, 

(pV(z)Pt(O))* 

<p;, + (z-)p;* + (o)>, 

a,b = l,2. 



-+(Y b Uz)Va-(o)) i 

=+<v; + (z)¥ , fl+ (o)>, 

-(P2L(z)PU(0)>,= 



— ) 

z 

_ Sqfe 

— — ? 

Z 

Z 

Safe 
_ ) 
Z 



(4.1) 



A point r = (n,^) <E R 2 in the two-dimensional Euclidean 
plane is hereby parametrized by the complex variables z := 
n + ;>2 and z '■= r\ — ir%, respectively. 

The disorder-average-smeared DOS in the MDH model is 
to be calculated by adding to the critical theory (I3.17> the per- 



turbation L 



£-ir| 



(e - CTl)o e _ iT1 where 



'21- 



'21- 



^124 
'^12- 



^21 + ) 

■^21-) 



(4.2a) 



and we have defined the composite operators 

^12+ := V1+V2-, ^2i+:=V*-V2*+. 
Fu- : = V2-V1+, 721-—V2+V1-, ( 4 - 2b ) 
in the fermionic sector [see Eq. (!3.15aH . and the composite 



operators 



>12H 



'12- 



P'l+P?- 
P2-P1+, 



'21- 



'21- 



PtP 2 *+, 



in the bosonic sector [see Eq. ( !3.15bH . Furthermore, in the 
definition of the DOS (I3.4ai the operator to be averaged over 
is now given by 



O v = +e 
+e 



-2*2 



'12+" 



'12- 



'21- 



'21- 



■^12H 
^12- 



^21H 
^21- 



(4.3) 
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Observe that O v has the same anomalous scaling dimension 
as the energy perturbation £ _ jT[ . 

Our strategy is, first, to derive a differential (Callan- 
Symanzik) equation obeyed by the DOS that encodes its 
change under the infinitesimal rescaling of the lattice constant 



(l+dl)a. 



(4.4) 



Second, we must solve the Callan-Symanzik equation obeyed 
by the DOS. 

In this spirit, we modify the critical action 5* defined by Eq. 
A3. 17t by adding to it the energy perturbation and any other 
operator o l needed to ensure invariance under infinitesimal 
rescaling of the lattice spacing. In other words, we seek a 
family of local composite operators {c\} such that 



S := S* - 



I 



d 2 r 



-2 



(4.5) 



is invariant under infinitesimal rescaling of the lattice spacing. 
The dimensionless coupling constants \ are scale (/) depen- 
dent and the anomalous scaling dimensions x v of the compos- 
ite operators O y are to be calculated at criticality. Invariance of 
the right- hand side under infinitesimal rescaling of the lattice 
spacing is equivalent to demanding that the family of com- 
posite operators labeled by e and I is closed under the OPE. 



Although the range of I is always infinite (for one, composite 
operators can always be point split and Taylor expanded), per- 
turbation of a unitary CFT simplifies greatly in that the fam- 
ily {o t } contains no more than a finite number of relevant 
(0 < x x < 2) operators. In practice, it is then sufficient to trun- 
cate the family {0^} to this finite subset of relevant operators 
if one is only interested in the effect of an arbitrary small per- 
turbation to criticality. For a nonunitary CFT, as is Eq. 13.17k 
the situation is dramatically different, however, since £ _ !T | 
induces under repeated OPE infinitely many relevant {x x < 2) 
operators. After selection of the family {o t }, the DOS v(e) 
obeys the Callan-Symanzik equation 



IT' 



P 



_ A 

^ dl 



(4.6) 



We then go on to solve Eq. ( 14. 6> in some approximate scheme. 



A. Operator content at criticality and RG equations 



At criticality and through OPE with itself, the energy per- 
turbation E _ (T1 induces local composite operators of the form 



o m , n = (^ 2+ p 2+ (^ 21+ p i+ (*i2 + r 2+ (« 21+ )" 2i+ (^,2-r 2 - (^-p- K 2 -p- K-p- 

x exp [2 (tn n+ +«?2i+ + "12+ + n 2i + — m u- — m 2i— — /2 i2— — w 2i— ) ^2] ) Vm,n G N 4 . 
It can be shown along the lines of Refs.l36landl37lthat O m n carries the anomalous scaling dimension 

x m,n = i }n \2+- m \2-) 2 + { m 2\+- m 2\-f + \ n V^ 



'12-1 



'21+- 



'21-1 



[('«12+ - '"12-) + ('"21+ - w 21-) + ("12+ - "12-) + ("21+ - «21-)] 2 , 



r 



(4.7a) 



(4.7b) 



which should be contrasted with its engineering dimension 

m i2+ + OT i2-+ m 2i+ + m 2i-+"i2+ + »i2-+»2i + + "2i-- Ob- 
serve the difference in the contributions to the anomalous scal- 
ing dimensions between the fermion and boson sectors. Fermi 
statistics demands point splitting when taking the product of 
Grassmann variables. This leads to the quadratic dependence 
on components of m in Eq. ( I4.7bi when g A = 0. Boson 
statistics accommodates arbitrary positive integer powers of 
complex-valued fields. This leads to the linear dependence on 
components of n in Eq. j4.7^ when g A = 0. The effect of 
disorder is thus much more potent in the boson sector than it 
is in the fermion sector since the positive term, which is linear 
in n when g A = 0, is always dominated by the negative term, 
which is quadratic in n when g A > 0, for sufficiently large \n\. 
The boson sector of the SUSY theory conspires with the un- 
boundedness of the exponential factor exp(2<t> 2 ) that enters in 



zero modes to produce infinitely many operators with negative 
anomalous scaling dimensions. 

Alternatively, the difference between fermion and boson 
sectors can be deduced by noting that the boson sector is re- 
sponsible for the normalization factor in correlation functions. 
Suppose that we perform an expansion of some correlation 
function in powers of exp(±2<J> 2 ) without using the replica or 
SUSY method. There are contributions to this expansion that 
originate from the numerator and from the denominator of the 
correlation function. Perform term-by-term disorder averag- 
ing. It was shown in the context of the 2D random phase XY 
model in Ref. |^] that (i) the counterpart to the contribution 
from the boson (fermion) sector to the anomalous scaling di- 
mension j4.7bt arises from disorder averaging over the contri- 
butions from the denominator (numerator) to the perturbative 
expansion in powers of exp(±2<J> 2 ) and (ii) disorder averaging 
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over the perturbative expansion in powers of exp(±2<J> 2 ) leads 
to the breakdown of perturbation theory. The lesson from the 
2D random phase XY model is that we cannot truncate the 
family {o m Tl } to some finite subset. We shall nevertheless 
argue that the family {Om. n } can be truncated to an infinite 
but countable subset when computing the DOS for the MDH 
model. 

To this end, we will again rely on the lessons to be drawn 
from a perturbative expansion of correlation functions in pow- 
ers of exp(±2<t> 2 ). To Mh order in powers of exp(±2<t> 2 ) the 
most singular contributions to the disorder average are associ- 
ated with terms proportional to exp(±2Mt> 2 ). For given 



labeled by < N £ N and ±. Be aware of the sign differ- 
ence on the second lines of Eqs. ( 14. 1 1 ai and d4.11b> , respec- 
tively. The choice {A±. N } C {o m n } is dictated by the alge- 
braic structure gl(2|2) of the critical theory ( I3.17t . Indeed, by 
construction, the left-hand side of 



N = m n+ +m 2l+ + n l2+ + n 2l+ 
-m l2 _-m 2l _-n l2 _-n 2l _ 



(4.8) 



the most relevant operators of the form (I4.7a> are given by 



N-n e ± 2 N® 2 



(®12±J [*21±, 

N > n > 0, 
N - 1 > n > 0, 

(^21±j («12±J ( S 21±j e 2 : 

N—l>n>0, 

{?n ± )(?2 i± )M(*2 i± ) N - 2 -"e ±2N ^, 

N - 2 > n > 0, 



(4.9) 



„±2A"I> 



' «12 ± + «21 ± T^12 ± T^21 ± ) Ar 



„±2/V<I> 



N\ 



«12± + «21±) 



e ±2N<£ 2 



^ N ^T~ (^12± + J 2 1±) («12± + »21±) 



N-l 



e ±lm 2 N _ 2 
+ N(N-l)^-f l2± T 2l± (« 12± + « 21± ) 



+ ■ 



(4.12) 



with the (2 x 4A0-fold degenerate anomalous scaling dimen- 



x N =N-^N z 

71 



(4.10) 



The (2 x 4A0-fold degeneracy of operators with the anomalous 
scaling dimension ( 14. 10> suggests that we truncate the family 
,7a> to the subset 



-:A ! 



-;N 



h2N<t>^ 



[(a 



12^ 



-N(f l2+ + T 2l+ ){'B l2+ + *2i+F~ 1 
+N(N -\)j n+ T 2l + (V n+ + V n+ ) N - 2 } 



(4.11a) 



-2W*, 



[(S 12 _ + « 21 _)' 



+N(!F n _ + T 2l J){<B n _ + <B 2l _f 1 

+N(N-l)T l2 _T 2l _(* l 2- + 'B2i-) N - 2 ], 

(4.11b) 



is a SUSY singlet [though not a gl(2|2) singlet]. Contri- 
butions to the binomial expansion of the left-hand side are 
related to higher-dimensional irreducible representations of 
gl(2|2)i22i Similarly, by taking the Mh power of the opera- 
tor O v in Eq. ( I4.3l l we can construct the counterparts Ay±-jj 
to A±;N that obey the same closed OPE's and share the same 
anomalous scaling dimensions. The consistency of the choice 
( 14. 1H is guaranteed by the fact that the OPE's (see Appendix 
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a + . N {z,z)A + . N ,{w 1 w) 
A_. N (z,z)A_. N ,(w,w) 
J4 + . N+ (z,z)A_. N _(w,w) 

M N _{Z,Z)S4 +N (w,w) 



,-x -x , ( N + N' 
N 



,-x -x , I N + N' 

W \ X N+N' X N X N' ( ' 



A_. N+N ,{w,w) + - 



|c-vi.'| % --' v + x,v - % + 



AT_ — 1 



^-;AL-iV,( W > W ) + ' 



< JV,JV' G N, 



< N,N' e N, 



(4.13a) 



(4.13b) 



<N_ <N + EN, 



(4.13c) 



< N, < N_ G N, 



(4.13d) 



r 



are closed within the subset {Ji±.^} C {o,^^}. Equation 
d4 . 1 3 1 is the most important technical result of this section. 
With the choice ( 14. lit . Eq. j4.5t is approximated by 



S = s *-L E d2 rAsM s > 

s=±N, = l J 



~ x s;N< — X N, — "s 



(4.14) 



Furthermore, by neglecting processes of "charge annihilation" 
as encoded by Eqs. ( 14. 1 3cl > and J4. 1 3di and their extensions 
when N + = N_ in favor of "charge fusion" as encoded by Eqs. 
(14. 13ai and J4. 1 3bi we obtain the counterpart 



E E 



dY 



S ;N S 



dY„. 



dl 



y s=±N<=\ ^1 dY s . N 



2 — X s;N ) Y S . N 



v(e), 



N s -1 



(4.15) 



to the Callan-Symanzik equation (14. 6> . At last, initial condi- 
tions on the coupling constants (fugacities) Y yN are indepen- 
dent of s = ± and given by 



Y s , Ns (l=0) = V^ 5 W 

a x s:< 



(4.16) 



The justification for neglecting all annihilation processes is 
that this approximation preserves the freezing transition that 
characterizes both the multifractal spectrum of zero modes 38 
and, as we shall see below, the DOS. This approximation is 
expected to capture the leading singularity of the DOS close 
to the band center. Furthermore, this approximation is con- 
sistent with neglecting renormalization effects on g A induced 
by the energy perturbation (see Appendix|XJ. At this level of 
approximation, the distinction between positive and negative 



charges is irrelevant to the Callan-Symanzik equation obeyed 
by the DOS. Equation ( 14. 151 thus simplifies to 



o= E^av--*i v(e) 

\N=\ U1 N ) 

Py =(2-x N )Y N + n £ ( " )Y N ,Y ! 



N'=\ 



x N =N-^N\ 
% 

Y N (l = 0) = £ - T ^8 LN . 
a i 



(4.17) 



It can be shown that Eq. J4. 17b is also obeyed by the DOS 
for a single flavor of Dirac fermions coupled to a real-valued 
random vector potential once all annihilation processes in the 
RG analysis have been neglected. Furthermore, essentially 
the same P functions for the fugacities in Eq. fl4. 17l > encode a 
RG analysis of the probability distribution for the normaliza- 



tion 



d 2 re- 2 ^) 



(4.18) 



of the squared zero mode amplitude in the MDH model pro- 
vided one replaces the anomalous scaling dimension x N = 
N — (g A /n)N 2 by — (g A /n)N 2 . This is most easily seen af- 
ter (i) performing the decoupling transformation J3. 14l > on Eq. 
( I3.4at applied to H MDii , (ii) isolating in the source term all 
terms that couple to exp(— 2<J> 2 ), (iii) ignoring the RG con- 
tributions from the SUSY matter fields, and (iv) identifying 
the energy e with L~ z a where L is the system size and z A is 
the dynamical scaling exponent that relates scaling with re- 
spect to L to scaling with respect to £. We have thus explic- 
itly related the transformation law obeyed by the DOS under 
rescaling to the multifractal properties of zero modes. Simi- 
lar RG equations were derived by Carpentier and Le Doussal 
within a replicated approach in Ref s . l35l andl44l in the context 
of the 2D random XY model and the freezing transition of the 
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multifractal spectrum of zero modes in the random vector po- 
tential problem, respectively. We now turn to the analysis of 
the Callan-Symanzik equation (I4.17t . 



B. From the Callan-Symanzik equation to the KPP equation 

The Callan-Symanzik equation J4.17t reflects the transfor- 
mation law under infinitesimal rescaling of the ultraviolet cut- 
off (lattice spacing) of the SUSY partition function 

Z := /'©[a,a ) 4» l! 4» 2 ,V* ) V J P'*,P']e-^ V£ , 



-2 



A. 



(4.19) 



where L % is given by Eq. (13. 171 and A n is, say, defined by 
the right-hand side of Eq. ( 14.1 lbi . As noted earlier [see Eq. 

633 ], 



(4.20) 



to leading order in an expansion organized by the relevance of 
the operators occurring in OPE's. It is then tempting to think 
of the coupling constants (fugacities) Y n as being related to 
the cumulants of some probability distribution. 41 Initially, this 
probability distribution is very narrow since only the first cu- 
mulant is nonvanishing. However, as the lattice rescaling pa- 
rameter / is increased, the cumulants grow; i.e., the probability 
distribution becomes broader. In the infrared limit / — > °° the 
ratios Y n J (3^)" diverge. The asymptotic RG flow of the fugac- 
ities thus offers little information on the bulk of the underlying 
probability distribution. Hence, it would be desirable to recast 
the RG equations obeyed by the fugacities Y n in terms of their 
generating function and to follow the RG flow of the gener- 
ating function instead of that of the cumulants it generates so 
as to obtain some information on the bulk of the generating 
function for the fugacities. Equation ( 14.201 suggests that the 
generating function for the fugacities takes the form 



®(^l) = t^Y n (l), 

/; = 1 



with the initial condition 



(4.21a) 



(4.21b) 



by Eq. ( 14. 171 . Differentiation with respect to the lattice rescal- 
ing parameter / gives, with the help of Eq. J4. 17i . the second- 
order nonlinear partial differential equation 



3,0 



In turn, reparametrization of and Z, 



©£,/)=:- G(y,l)-l 
K L 



© + 7I0 2 . (4.21c) 



% =: -e-y, (4.22a) 



d,G= (a v + ^a; )g+2g(g-\), 



G(y,0) 



2 a x i~ 2 



(4.22b) 
(4.22c) 



The convection term on the right-hand side is induced by the 
RG flow of the matter fields in the bosonic sector. This term 
is absent from the RG analysis of the normalization d4 . 1 8 i of 
the zero mode. At the level of the generating function, the 
convection term can be disposed of by the change of variable 
consisting of a Galilean boost and of a rescaling 



G(x,t):=G{^/^x-t/2,t/2), 
after which 

d,G=^d 2 x G + G{G-l), 



GM) = i-^V? ex P 

1 a i 



! exp 



(4.23a) 



(4.23b) 




(4.23c) 



The error caused by exponentiation in the last step should 
be comparable to the one caused by neglecting the renormal- 
ization of g A induced by the energy perturbation. Equation 
(I4.23bl is known in the mathematical literature as the KPP 
equation with the forcing term 



F(G):=G(G-l) 



(4.24) 



when* G K and < G< 

Trading 0(^,/) for G(y,l) [or G(x,t)] offers the following 
advantage. It is possible to interpret G(y,l) as the generating 
function of a random variable Z with the (scale-dependent) 
probability distribution P(Z,l) through the definition 

G(y,l) =: (exp(-e- y Z)) p(zl) (4.25a) 

with the initial condition 

P{Z, 0) = 8(Z - %Y l (0) /2) . (4.25b) 

This definition is consistent since the probability distribution 
P{Z,l) is normalized to 1 initially and since its normalization 
is invariant under the RG flow. In turn, we can reinterpret the 
fugacities Y„ as the moments of the (scale-dependent) proba- 
bility distribution P(Z,l) with the initial conditions 



(4.25c) 
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Equation (I4.25c> should be compared with Eq. J4. 17t . 
This immediately leads to the identification of the (scale- 
dependent) random variable Z with the scale-dependent ran- 
dom energy (;t/2)[e(/) - ir\(l)]/a 



x, —2 



i.e. 



G(y,l) = 1 + -©£,/) 



exp 



_ y ne(l) - 
2 



-2 



P(Z,I) 



and, more importantly, 



Z(0)/ 



P(Z,l) 



e(0)/ 



p(z,/) 



-lyp 



(0 



typ 



(0) 



:= exp 



zm \ 



/•(z,/) 



(4.26a) 

) 

(4.26b) 
(4.26c) 



As we shall see in the next two subsections, it is possible to 
extract from Eq. ( I4.26ct a scaling exponent z A called the dy- 
namical scaling exponent that encodes the scaling of the typi- 
cal energy £ typ (/) in the infrared limit / — ► °° through 



z A := lim -In —, ; . 
A I e tyP (0) 



(4.27) 



Knowledge of the dynamical scaling exponent z A suffices to 
compute the DOS in the MDH model. 



C. Analysis of the KPP equation 



In the limit t — > °°, the KPP equation ( I4.23bl > possesses a 
solution in the form of a traveling wave w A (x — v A t) that is 
unique up to translation and whose velocity v A is solely con- 
trolled by the exponential decay of the initial data for large x, 
provided the initial data obey the conditions 



G(x,0) = 



1 , for x 



0, for x — > — °°, 



(4.28) 



which they do according to Eq. ( I4.23c^ More precisely, if we 
choose the origin of x such that 



limll-GM)]^-^ 75 *, 



(4.29) 



then 



LJZa 

2 V 71 



V2, 



and 



lim[l-w A (*)] 



lf A , for^<2, 



for & > 2, 

rV^ 1 , forv A >V2, 
ce'^, forv 4 =V2. 



(4.30) 



(4.31) 



It is a remarkable property of the KPP equation^i^ that the 
form of the forcing term F(G) in Eq. does not affect 

the velocity of the traveling wave as long as 

F(0)=F(1)=0, F(G)<0, 

F'(0) = -1, F'(G)>-1, 0<G<1. 

(4.32) 

As noted in Refs. l35l and l44l had we chosen a different RG 
scheme such as smoothened the hard-core condition in real 
space used to obtain Eq. ( 14. 151 . say, we would have obtained 
another forcing term in the KPP equation. Universality, i.e., 
the independence of the choice of RG scheme, is thus equiv- 
alent to finding a forcing term that obeys Eq. ( 14.321 from the 
point of view of the KPP equation. 

Equations J4.30l > and ( 14.3 11 allow us to perform the sum- 
mation over all fugacities in the generating function ( I4.23a> 
when / — > °°: 



1-G(*,2/) = 



-Ui) 



r VgIfr(x-2v A ') forv >y /2 



,-v5(i-2v A 



for v A = y/2, 



(4.33) 



as x — > oo. Although individual fugacities are relevant and 
grow without bounds in the infrared limit, a finite limiting 
value of the generating function can nevertheless be extracted 
if x scales linearly with / as a result of the traveling-wave 
asymptotic of G(x,2l). 



D. KPP equation and DOS 



We are ready to make an estimate for Eq. ( I4.26bi that justi- 
fies the definition of the dynamical exponent (I4.27> . To this 
end, we take advantage of the facts that, by Eqs. (I4.23a> . 
( !4.23cL and ( 14.331 . the generating function G(x,2l) can be 
thought of as a wave front which, initially, was chosen to 
coincide with the origin x = and is located at x — 2v A l in 
the long-time limit t = 21 —> °°, In the same way, G(y,l) = 
G(^%/g A (y + /),2/) is a wave front located at y = and 
y = (2y / g A /nv A — 1)1 initially and in the infrared limit / — > oo, 
respectively. This gives the estimate 



In 



Z(0)/ 



P(Z,i) 




G(y,0)-G(y,l] 



G(x,Q) 



-G{x + ly/n/gA,2l) 




0(lnZ), 



(4.34) 
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where the logarithmic corrections are universal (see Ref. |5 
It is then tempting to identify 



£typ(0 

MO) 
whereby 



exp < In 



z(/n 

Z(0)/ 



oc e A 



(4.35) 



d\m t 
hm — — ^ 

(-k*> dl 



iorg A < 2k, 



(4.36) 



defines the quenched average value z A of the dynamical expo- 
nent, as a coupling constant that enters the Callan-Symanzik 
equation obeyed by the DOS, which we here assume to be 
self-averaging and to obey some sort of scaling ansatz. This 
conjecture is bolstered by the fact (see Refs.l45landl46T) that 
the distribution of lnZ(7) is self-averaging in the infrared limit 
/ — > oo as we had hoped for. Consequently, it is plausible that 
£ typ is the only relevant coupling constant that enters 



typ 



the Callan-Symanzik equation for the (self-averaging) DOS 



d_ 



v(e). 



(4.37) 



forg 



< 2k.) A rigorous 
J4.37I and ( 14.171 are identical 
-> oo at some fixed energy e is 



(Observe that 2 — z A = l—g A /n- 
proof that the solutions to Eqs. 
in the thermodynamic limit L 
beyond the scope of this paper. 

The mechanism by which solutions to the KPP equations 
are traveling waves causes the dynamical exponent z A to un- 
dergo a freezing transition when the disorder strength g A 
reaches the critical value 2%. The freezing transition of the 
dynamical exponent reflects itself by a nonanalyticity of the 
multifractal spectrum of zero modes as was shown in Refs.l38l 
and|4(J It also reveals itself through a nonanalytic dependence 
on g A in the power-law behavior of the DOS since integration 
of Eq. J4.37i yields 



v(e) 



v(ej 



for?- < 2k, 



(4.38) 



The energy scale e r . n is the one for which errors caused by ne- 
glecting renormalization effects on g A are of order 1. For 
weak disorder (g A < 2k), quenched and annealed average val- 
ues of the dynamical exponent agree since 2 — z A =x v For 
strong disorder (g A > 2k), quenched and annealed average 
values of the dynamical exponent differ since 2 — z A 



Bounds on the dynamical exponent that are consistent with 
Eq. ( 14.361 were established in Ref. ^3 The same power-law 
dependence for the DOS of a single species of Dirac fermions 
subjected to a real random vector potential was derived in Ref. 
'9 for weak disorder. In this simpler context, the transition 
from weak to strong disorder for the DOS was established in 
Ref. |^| using a replica variational method. We now turn our 
attention to the computation of the DOS in the HWK model 
for which a field-theoretical approach is lacking. 



V. DOS FOR /W 

The HWK model differs from the MDH model by the pres- 
ence of a complex-valued random mass M_(r) = [M + (r)}*. 
It was pointed out in Ref. |22| that, although g A is exactly 
marginal in the MDH model (13. 171 . the presence of a random 
mass renders g A marginally relevant (in the sense that it grows 
logarithmically with length) whereas g M remains marginal in 
a one-loop RG analysis. Guruswamy et al. computed the ex- 
act (3 functions for g A and g M , thereby confirming the one- 
loop RG analysis: 2 - Furthermore, they found the anomalous 
scaling dimension 



*e = C[l + 0(ft*)]-^ 



(5.1) 



for the energy perturbation in Eq. ( 13 . 1 5i by taking advantage 
of enhanced symmetries of the HWK model at the band cen- 
ter. Remarkably, the dependence on g A in Eq. ( 15. 1> is exact 
while that on g M can be obtained by perturbation theory. Fi- 
nally, they wrote down a Callan-Symanzik equation for the 
DOS of the form ( 14. 6t with a finite number (3) of coupling 
constants. The three coupling constants are e, g A , and g M . 
The corresponding p functions are 



p £ = (2-x E )£, p ?A = 



2k 2 



S M 



(5.2) 



respectively. Additionally, the anomalous scaling dimension 
x £ entering the Callan-Symanzik equation is given in Eq. ( 15. Q . 
Solving this Callan-Symanzik equation gives the diverging 
DOS dl . Il l in the vicinity of the band center. 

One can tune the HWK model to the MDH model by 
switching off the random mass perturbation, i.e., by taking the 
limit g M = 0. In this limit, the Callan-Symanzik equation used 
by Guruswamy et al. only agrees with the Callan-Symanzik 
equation J4.37l > for the DOS when g A < 2k, i.e., below the 
freezing transition. As we have seen in Sec. IIVI this discrep- 
ancy is caused by a failure to account for all the relevant cou- 
pling constants associated with powers of the energy perturba- 
tion. In this section, we are going to show that infinitely many 
relevant operators were overlooked by Guruswamy et al. in 
their RG analysis of Eq. d3 . 1 51 at the band center. We will then 
derive a generalized KPP equation for the generating function 
of these coupling constants. This generalized KPP equation 
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differs from the KPP equation J4.22i by a scale-dependent 
diffusion "constant" driven by the scale dependence of g A . 
Within an adiabatic approximation, we shall extract from this 
generalized KPP equation a Callan-Symanzik equation for the 
DOS of the form (I4.37> except for the presence of the run- 
ning coupling constant g A whose p function is given by Eq. 
\5.2\ . However small the bare value of g A , it always flows to 
infinitely strong disorder. The freezing transition must then 
necessarily affect the DOS in the HWK model. The solu- 
tion for the DOS in the HWK model is in fact given by Eq. 
dl.lt with the exponent K = 2/3 instead of K = 1/2. It is 
through the larger exponent K = 2/3 that the freezing tran- 
sition reveals itself as was noted by Motrunich et al*& A more 
detailed analysis of the Callan-Symanzik equation shows that 
the asymptotic regime JO with K = 2/3 can be preceded 
by crossover regimes characterized by the Gade singularity or 
by the power-law behavior of the MDH model provided bare 
coupling constants are fine-tuned. 

A. Nearly conformal field theory 

The chiral decoupling transformation j3.14t on the HWK 
model in the SUSY representation ( 13 . 1 3i is not as useful when 
g M > as it is when g M = since rotated SUSY matter fields 
remain coupled to the disorder through the random mass. Our 
starting point in this section will thus be the effective theory 
obtained from the HWK model at the band center after inte- 
gration over the Gaussian-distributed random vector and mass 
potentials 

Z:=/«[ ¥ *, ¥ ,r,P]^ /421 ( £0+ " OA+ ^ OM l 

2 

:= E (¥a-3v fl -+¥fl + 9v fl + 

a=l 

2 

°a=- L (v:-v„-+p;,-P„-)(Vi + Vi + +Pfo + p fo+ ), 

a,b=\ 
2 

M =- L (V«+V«-+P«+Pa-) (vt-Vb++Vb-Vb+) • 

a,b=\ 

(5.3) 

Here d := and d := jd_. Observe that the effective the- 
ory (15.31 is invariant under any permutation of the flavor in- 



dices. Flavor permutation symmetry is broken by the energy 
perturbation ( 13. 13cl >: i.e., the SUSY symmetry group in the 
presence of an energy perturbation is enlarged to the group 
GL(2|2) at the band center. Since we are after the behavior 
of the DOS in the immediate vicinity of the band center, it is 
essential to account for the enlarged SUSY symmetry group 
GL(2|2) at the band center. 

Fermionic interacting field theories can be considerably 
simplified by the method of bosonization. This is also the 
case of the SUSY representation J5.3i of the HWK model at 
the band center as shown by Guruswamy et a/.^i and we now 
summarize their results. We thus adopt Abelian bosonization 
rules that are summarized in TablesH HIVl S Abelian bosoniza- 
tion in the fermionic sector is performed by trading left, \|/*_, 
say, and right, Vf* + , say, moving fermions for the exponential 
of holomorphic cp„_ and antiholomorphic (p a+ scalar fields, 
respectively. Abelian bosonization in the bosonic sector was 
introduced in string theory^ii^ It is performed by trading left, 
P*_, say, and right, p* + , say, moving bosons for the expo- 
nential of holomorphic (p',_ and antiholomorphic (p' (+ scalar 
fields, respectively. Additional fermionic fields %* and % a 
must also be introduced when bosonizing bosonic spinors. A 
simplifying feature of the two interactions induced by the dis- 
order in Eq. \53\ is that they can be written as bilinears in the 
gl(2|2) currents that generate the GL(2|2) symmetry group. 
These currents are of two kinds. The first kind is made of 
bosonic currents denoted / ■ , J J ' b , and J 7 which are de- 
fined in tables U and HI] The second kind is made of fermionic 
currents denoted G ab+ , G ab+ , G ub _, and G ub _, which are 
also defined in Tables [Q and [H] Currents of the bosonic 
kind permute the flavors of right and left movers but do not 
change their statistics. Currents of the fermionic kind ex- 
change fermionic spinors of one flavor for bosonic spinors of 
another flavor, and conversely. The importance of the gl(Af |A0, 
NeN algebra for Dirac fermions coupled to random fields 
was emphasized in Refs. |3(J |^ and |3 The central role 
of the symmetry group GL(Af|A0 also appears in the quan- 
tum field theory description of the multivariable Alexander 
polynomial^ 

After Abelian bosonization, the holomorphic and antiholo- 
morphic scalar fields (p a± and cp', ± enter Eq. (15.31 through the 
combinations 

§a '= <P«-+<P a +, ¥a ■= + a=l,2 

(5.4) 

only (see TablesllTlland lIV> . 
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TABLE I: Abelian bosonization rules for gl(2|2) matter fields and currents: flavor-diagonal components of the currents. The scalar fields (|) a , 
(p fl - , and (p n+ are related by (|) fl (z, z) =<p a - (z) + <Pa+ (z) ■ The scalar fields §' a , (p',_ , and (p' I+ are related by §' a (z,z) = cp„_ (z) + <p' a+ (z) ■ The 
lattice spacing a is used as a short-distance cutoff. The magnitude of the complex number u is arbitrary as it reflects the nonuniversality of the 
short-distance regularization. The phase of the complex number u is fixed by the boundary condition obeyed by the scalar fields at infinity, 
which, in turn, is uniquely fixed by the SUSY of the theory. Demanding that the scalar fields vanish at infinity amounts to fixing the phase of 
u to be — ir/2. We will choose u = —i for convenience. 

SUSY matter Bosonized matter SUSY currents Bosonized currents 



¥„- 


u +l/2 a -l/2 e -«p„- 


ha ■= ¥„-¥„- 




(+;'3(|> a ) 




¥„- 


u -l/2 a -l/2 e+i% . 


ha := ¥«+¥«+ 




(— l'5«t>a) 




¥„+ 


( ; ,*) + 1 /2 a - 1 /2 e +/<P n+ 


ha • P(7— Pfl— 








¥ fl+ 




ha := Pfl+Po+ 








K- 


„-i/2 a +i/2 e -i<_ (i -a Xa) 


G aa+ '■= Pfl-¥o- 


II" 1 


e -'(<pi--<p<,-) 




P a - 


„+i/2 a +l l 2 e+<-{id%i) 


G aa+ '■= Pfl+¥o+ 




i e +'(f>;,+-<fw) 




p:+ 


{u*)- x l 2 a +l l 2 e +i ^{idx a ) 


G aa- '■= Pn-¥a- 


11 


e +'(<p»--<p«,-) 




p,+ 


{u*)+ l l 2 a+ l l 2 e- i ^{idx* a ) 


G aa- '■= Pa+¥a+ 


u* 




{>h* a 



TABLE II: Abelian bosonization rules for gl(2|2) matter fields and currents: flavor-off-diagonal components of the currents, i.e., a^b. 
The scalar fields § a , (p a _, and (p fl+ are related by <]) a (z,z) = <Pa-(z) + 9a+(z). The scalar fields (])',, (pj,_, and q>' a+ are related by 0a(z,z) = 
9a_(z) + 9a+(z). The lattice spacing a is used as a short-distance cutoff. The magnitude of the complex number u is arbitrary as it reflects 
the nonuniversality of the short-distance regularization. The phase of the complex number u is fixed by the boundary condition obeyed by the 
scalar fields at infinity, which, in turn, is uniquely fixed by the SUSY of the theory. Demanding that the scalar fields vanish at infinity amounts 
to fixing the phase of u to be — jt/2. We will choose u = —i for convenience. 

SUSY matter Bosonized matter SUSY currents Bosonized currents 





„+l/2 a -i/2 e -i%- 


ha 


= Va-Vt- 


Q -i e +/(<p„_- 96 _) 


¥<,- 


u -l/2 a -l/2 e +i%- 


h a 


= ¥ fl+ ¥fc+ 


Q -l e -'(<(>«+-96+) 


¥„+ 


( M *)+V2 a -V2 e +i<P a+ 


J ba 


= P fl -PL 




¥«+ 


( M *)-V2 a -l/2 e -i<P« + 


J ba 


= P a+ P^ + 


ae-^-^Hiha^hb 


PS- 


u -l/2 a +l/2 e -,<_( ( -9 z J 


G ba 


f := P«-¥ fc - 


„-l e -^-^)(0 Xfl 


p fl - 


M +l/2 Q +l/2 e +/<-(/3 x *) 


G ba 


f := K+W b+ 


(M * r i e +<W + -<P 4+ ) (! -3 Xa 






G ba 


- ■= P«-¥/*- 




p a+ 


(u*) +1/2 a +1 / 2 e-^ + (ih* a ) 


G ba 


- := P a+ ¥fe+ 





Z= / ©[^,(j)^x:,X«]exp( - / i-^L 



2n 



(5.5a) 



£, = 1 



a=\ 



2 

L 



l + iK e -i(*a-Va) 



(3xS)(3Xay 



8m 



ajtb 



(5.5b) 



(5.5c) 



Here, we have introduced the "center of mass" E + and the "relative coordinate" E_ for the scalar fields 

2 2 

E + :=E + E', E_:=E-E', £:= J> fl> := JX- 

a= 1 a=\ 



(5.6) 
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TABLE III: Abelian bosonization rules for mass and vector potential current-current interactions: flavor-diagonal contributions. Columns 1 
and 2 are the contributions to the mass perturbation and columns 3 and 4 are the contributions to the vector perturbation. With the choice 



—i in Tables HlandlrIL 



I, 



SUSY matter 



: ll[U* 



-1. 



Bosonized matter 



SUSY matter 



Bosonized matter 



J aa' J aa 
+ G aa+ G a 

—G„„_G„ 



+(3<W(3W 

» e -^-«(;a X£ ,)(/3 x *) 



-u(u* 



u aa u aa 

+j f 

1 aa aa 

+/ / 



-(3* B )(3*«) 

+(^J(3C) 
+W a )(34»«) 



TABLE IV: Abelian bosonization rules for mass and vector potential current-current interactions: flavor-off-diagonal contributions, i.e., a^b. 
The lattice spacing a is used as a short-distance cutoff. Columns 1 and 2 are the contributions to the mass perturbation, and columns 3 and 4 



are the contributions to the vector perturbation. With the choice u = —i in TablesHlandlrIL 



-1. 



SUSY matter 



Bosonized matter 



SUSY matter 



Bosonized matter 



+ J ba J ab 
J ba J ab_ 

+ G ab+ G ab- 



-G 



ba- 



ba+ 



+a -2 e +i(*.-fc) 

-a+ 2 e+W°-KXdxt)(dz b )(ht)(h a ) 
+i l - 1 i,*e- i ^-^(idx b ){idx* b ) 
-u(u*)-ie+W«-M(®X*)(idx a ) 



^aJbb 

-/' V 
J aa J bb 

^■*aJbb 



-(3<M(§<fe) 

-Wa)(Wb) 



After Abelian bosonization, the decoupling transformation P.14t takes the very simple and useful form 



(5.7) 



As pointed out by Guruswamy et a/.^lEq. (15. 7> is crucial 
to establish that the HWK model is, in essence, a nearly con- 
formal field theory at the band center. Equation d5 .71 follows 
from the peculiarity of the Lagrangian L by which the center 
of mass E + only enters (in L^) through the linear term 



(5.8) 



To put it differently, the term (15.81 is the only one in L affected 
by the chiral transformation 



a = 1,2. 



(5.9) 



All other terms in L and, in particular, all terms in L 2 are left 
unchanged by the chiral transformation ( 15. 9> as is best seen 
after trading in L the four coordinates fl and (a = 1,2) 
for E + , E_ and the relative coordinates a a+ := O a + a', and 
& a _ := a„ - & a whereby 4> fl =: (1/2) + o fl and =: (E'/2) + 
O^,. The role of the center of mass E + is thus to enforce a 
global constraint on the coordinates a and §' a , a = 1,2 of the 
gl(2|2)-graded Lie algebra. 

We are mostly (if not only) concerned with the evaluation of 
the expectation value of the product of vertex operators (up to 
Klein factors that do not matter for the calculation of anoma- 



lous scaling dimensions) 



(id Xl ,), (5.10) 



where I and l' are collective labels for flavor and space coordi- 
nates, should be understood as any of the four options 

(idx x ,), (ify,), (rd%*,), and (id%*,), and a p a[, € R. Under the 
chiral (additive) transformation law ( 15. 9> the transformation 
law of a product of vertex operators is multiplicative 

X^[<LCXnXj- (5.H) 

The usefulness of Eqs. ( I5.7i . ( 15. 9> . and ( 15.1 It is the following. 
If the expectation value of the left-hand side of Eq. ( 15.1 1> is to 
be evaluated with the partition function J5.3i at given g M and 
g A , then the expectation value of the right-hand side of Eq. 
( 15. lit is to be evaluated at the same g M but different g A = 0. 
Let us denote by ((■••)) . , the expectation value with the 
partition function J5.3I . We have just proved that 



(5.12) 

Let us denote by g [§ a Aai%*»%a\ me right-hand side of Eq. 
(I5.1H . We also write ((■••)} „ for the expectation value 
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with L in the partition function ( 15. 5> replaced by L x with 
g A = 0. We then define the perturbative expansion of the 
expectation value (g [^ a Aa,%t,X a })x,,g M ,0 around the "unper- 
turbed theory" L , with g A = by 



ISM' 



r (n) 



d 2 r Y 
~2~K 



d 1 r„ 



-<£[<!>« X%a,X J 



2k 

x L 2 (ri)---L 2 {r n ) 



t v g M ,0 

(5.13) 



To go from Eq. (15. 121 to Eq. ( 15 . 1 3i we made use of the fact 
that the perturbation L 2 does not couple to the center of mass 
E + and is thus invariant under the chiral transformation d5.9i . 
Nearly conformal invariance is a property of the theory de- 
fined by Lagrangian L Y with g A = as L i then reduces to 



two independent copies of the N =1 limiting theory for which 
nearly conformal invariance was demonstrated in Ref.|3 Be- 
fore we characterize in more detail what is meant by nearly 
conformal invariance, it suffices here to say that any perturba- 
tive contribution to the right-hand side of Eq. (15. 12i is to be 
evaluated with the following rules, (i) All fields obey Wick's 
theorem, (ii) Any two-point function between £_ and the rel- 
ative coordinates entering L 2 vanishes. If so, perturbation the- 
ory in powers of L 2 does not renormalize the dependence on 
g A of the expectation value of the right-hand side of Eq. (15 . 1 3i 
obtained to zeroth order in powers of L 2 . We have thus proved 
that the exact dependence on g A of any functional of the form 
(15. lOt is already obtained to zeroth order in perturbation the- 
ory in powers of L 2 around the unperturbed Lagrangian £j 
with g A = once the chiral transformation ( 15. Ill has been 
performed. 

Under nearly conformal invariance of the theory defined by 
the Lagrangian £ j with g A = 0, we understand the following 
properties demonstrated in Ref.l27l 



(i) Nonvanishing two-point functions between the fields ty a , §' a [see Eq. 15.41 1. their duals 

<t>„ : = 9a- - %+, := 9a- - <P«+> 



C[-8 < *-(&iC/*)]ln 
^[+8 ah -(g A ^/n)]\n 

-(*aC7*)i° 



fe M C/(27T)] 2 C 2 ln 2 
fe M C/(2^)] 2 C 2 ln 2 
-fe M C/(27T)] 2 C 2 ln 2 



and the auxiliary fermions %* and % a are 53 

(♦ a (z,f)Mo,o)) Xiftf)gA 

(Za(z,z)Zi,(0,0)) £iigMi0 = -(x a (z,z)xt(0,0)) L[ , gM , = -8 ab ^\n 



zz 



= C [+8 ab + (g^/n)] In ( ^ ) + fe M C/(2^)] 2 C 2 ln 2 ( g 



a 

c)] 2 C 2 ln 2 (§ 



(5.14) 



(5.15a) 



(5.15b) 
(5.15c) 



(a = 1,2), respectively. There are nonvanishing correlation functions between (])'s and their dual (j)'s but we will not need them. 
Here, the ellipses means that short-distance contributions such as 8 functions have been omitted. Observe that Eqs. ( 15. 15al and 
d5 . 1 5bl > apply not only to L l with g A = but also to L with g A > 0. However, we stress that Eqs. ( 15. 15al and ( 15. 15bi were first 
computed with g A = 0. The full g A dependence was then uncovered through the chiral transformation d5.9l >. Another important 
point to keep in mind is that the squared logarithmic term enters with a negative sign in the two-point functions for the fields a 
and (j)',, while it enters with a positive sign for the dual fields <j) a and §' a . 

(ii) Many-point correlation functions among (j) a , §' a , (j) fl , §' a , %*, % a can be reduced to the two-point correlation functions d5 . 1 5I > 
with the help of Wick's theorem. Applicability of Wick's theorem is specific to L l with g A = 0. 



The (3 functions in Eqs. d5.2i for gA and gM follow from Eqs. d5. 15al > or d5. 15bi by demanding that these two-point 
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TABLE V: Abelian bosonization rules for the eight Cooper terms 
entering the energy perturbation 14.21 . The fields § a , (p fl -, and (p fl + 
are related by § a {z,i) = 9a- (z) + 9a+(z)- The fields §' a , y' a _, and 
<f/ a+ are related by <S?' a (z,z) = 9a-(z) +9 a +(z). The fields <j) a , 
and (p a+ are related by <j> a (z,z) = 9a-(z) — 9a+(z)- The fields 
(p„_, and (pj, + are related by <|>q(z,z) = cpj,_ (z) — 9 a +(z). The lattice 
spacing a is used as a short-distance cutoff. 



SUSY 



Bosonized Cooper terms 



9T-95+ 
Yl+9l- 

Pl+PI- 

V 2+ ¥i- 

¥ 2 -Vi+ 

P 2 +Pi- 

P 2 -Pi+ 



a 1 exp 




} 
{ 



functions be invariant under rescaling 04. 4> of the lattice spac- 
ing a. It is the presence of the squared logarithm terms in 
Eqs. (15 . 1 5i that spoils the conformal symmetry of the MDH 
model, i.e., the theory with gA > and gM = 0. The damages 
caused to conformal invariance are minimal, however, since 
it remains possible to use Wick's theorem perturbatively and 
since the dependence on gA in correlation functions is exact to 
lowest order in perturbation theory. 



B. Operator content at near criticality 



We have seen in Sec. lIV Al that the energy perturbation ( 14.21 
induces local composite operators given by Eq. ( I4.7ai with the 
anomalous scaling dimensions (I4.7bl in the MDH model. We 
are going to show that the dependence on g M of the anoma- 
lous scaling dimensions of these operators can be calculated 
perturbatively in the HWK model whereas the dependence on 
g A is, up to the positive multiplicative factor i^ 2 , exact in Eq. 
( I4.7bl . Hence, it is imperative to account for the most rele- 
vant of these operators in a consistent way when writing the 
Callan-Symanzik equation for the DOS. 

To begin with, we must bosonize the eight Cooper terms 
entering in the energy perturbation d4.2i . This is done with 



J 



the help of Table [V] from which we infer that only the four 
linear combinations 



and 



0_ 



(9i + 9 2 )±(9i-9 2 ) 



(4>i + 4>2)± (91-92) 



(5.16) 



(5.17) 



enter in the vertex operators representing the Cooper terms. 
Observe that the energy perturbation depends on the dual 
fields after Abelian bosonization. 

Permutation symmetry of the flavor indices at the band cen- 
ter simplifies considerably the evaluation of two-point func- 
tions between 0's. In particular, two-point functions between 
scalar fields and their dual cancel each other. Alternatively, 
this is so because the conformal spin of the energy perturba- 
tion vanishes; i.e., there are as many left as right movers in any 
of the monomials making up the energy perturbing operator. 

With the help of Eqs. ( 15. 1 5al > and ( 15 . 1 5bl > we can evalu- 
ate exactly all two-point functions listed in the first column of 
Table IVl! For the purpose of computing scaling dimensions, 
these exact two-point functions are only the zeroth-order con- 
tribution to the expectation value of products of vertex oper- 
ators in perturbation theory in powers of X 2 . They are given 
in the second column of Table IV1I Observe that the depen- 
dence on the squared logarithm and hence on g A (through 
the chiral transformation) arises solely from the dual sector 
of the theory as the contributions from the nondual sector al- 
ways cancel. This is consistent with the fact that the random 
vector potential only couples to the dual field. In the present 
formalism, this follows from the fact that the center of mass 
drops out in the difference of any of the coordinates (j) f , or 9',. 
The dependence on g A extracted from Table IVII is exact for 
the scaling dimensions of vertex operators in that perturbation 
theory only changes the linear dependence on the logarithm 
through corrections that depend solely on g M . Observe that 
vertex operators corresponding to products of (3's follow from 
two-point functions with prime fields. Hence, the contribution 
to the term linear in the logarithm that does not depend on g A 
has the "wrong sign." In the absence of a random mass, this 
wrong sign is corrected by accounting for the contribution to 
the correlation functions from the % sector [see Eq. d4.7bl and 
Appendix 151. In the presence of a random mass, the correc- 
tion to this "wrong sign" from the % sector is only partial [see 
Eq. d5~18t l. 



The central result of this section is the evaluation of the anomalous scaling dimensions x m n for the local composite operators 
defined by Eq. ( I4.7a> . To zeroth order in the perturbation theory defined in Sec. IV Al they are given by 
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TABLE VI: Evaluation of all two-point functions governing scaling dimensions appearing in OPE's of energy perturbation with itself. The 
ellipses refer to 8 function contributions to the two-point functions. Although these two-point functions are exact, they yield an approximation 
to the calculation of scaling dimensions. Nevertheless, the dependence on g A of scaling dimensions is exact. The dependence of linear 
logarithmic term on g M that does not depend on g A is modified by residual interactions in L when computing scaling dimensions. 
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as the expectation value of any product of vertex operators 
can be replaced by the exponential of the properly weighted 
sum of all 12 nonvanishing two-point functions listed in Table 
IVII Anomalous scaling dimensions are then extracted from 
the coefficient to the term linear in ln(fz/a 2 ). One should 
not forget the contribution to the anomalous scaling dimen- 
sion coming from the which is responsible for the third 
line in Eq. (I5.18> . What makes Eq. J5 . 1 81 remarkable is 
that even if O m n is not relevant at g M = 0, switching on a 
mass perturbation renders O m n relevant at sufficiently large 
length scales since g A is marginally relevant (in the sense that 
it grows logarithmically with length) whereas g M is exactly 
marginal. Moreover, higher-order corrections from perturba- 
tion theory do not invalidate this argument. The strategy to 
deal with this difficulty was described in Sec. II VI We are go- 
ing to generalize the RG analysis of Sec. II VI to account for 
the relevance of g A . As before we select the operators ( 14. 9> 
since they are the most relevant operators for given N in Eq. 
d4.8t although it must be noted that their scaling dimensions 
are not anymore degenerate for finite g M . The lifting of the 
degeneracy J4. lOi is negligible, however, since it is of order 
8m/ 8a wnere by g M does not flow but g A flows to infinity in 
the infrared limit. 



C. Callan-Symanzik equation for DOS 

The Callan-Symanzik equation obeyed by the DOS in the 
HWK model is given by 




Y NdY N ^M d g M ^S AdgA 



where the p functions for the fugacities are 

N-l f N 

V Yn = (2-x n )Y n +tz £ 

N'=l 



h )v(e), 

(5.19a) 



x N = Z ) [l + o(g M )}N- 



N' 



N 



(5.19b) 



the P function for the variance of the random mass vanishes, 
and the p function for the variance of the random vector po- 
tential is 



(8mQ 2 
2k 2 



(5.19c) 



[Remember that p and p , follow from demanding that the 
two-point functions Eqs. 15 . 1 5 al or J5 . 1 5bb be invariant un- 
der rescaling J4.4I of the lattice spacing a.] The function £ 
is given in Eq. ( 15.2k We made the same approximations to 
reach Eq. d5 . 1 91 as we did to obtain Eq. ( I4.17l i; i.e., we ig- 
nored renormalization effects induced on g M and g A by the 
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fugacities and we ignored annihilation processes when deriv- 
ing the RG equations obeyed by the fugacities. Furthermore, 
we work to lowest order in g M as g M does not flow whereas 
g A flows to infinity in the infrared limit. The infinite set of 
RG equations d5 . 1 9bt obeyed by the fugacities can be recast 
as a second-order nonlinear partial differential equation of the 
KPP type 
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The notable difference with Eq. j4.22> is the scale dependence 
of the diffusion term g A (l)C, 2 /n. 

We tentatively define the dynamical exponent z A for the 
HWK model by 



z A ■ = 



4 ^-C-ofe M ), for^C 2 >2:t. 



(5.21a) 

The dynamical exponent z A in Eq. J5 . 2 1 al would be the veloc- 
ity of the wave front of the Callan-Symanzik equation (15.201 
if the diffusion constant was scale independent. However, the 
dependence of the diffusion constant oc gj^{l) on the rescal- 
ing parameter / is linear whereas that of G(y,l) is exponential. 
We can thus make an adiabatic approximation by which the 
Callan-Symanzik equation for the DOS (I4.37> becomes 

Za<^-+pY^--(2-z a ) 



0: 



3e 8 *dg A 



v(e). 



(5.21b) 



The Callan-Symanzik equation d5 .2 II is the central result of 
our paper. It differs from the Callan-Symanzik equation of 
Guruswamy et al. in Ref.l27l through the freezing transition of 
the dynamical exponent in Eq. 15 . 2 1 ai . 



D. DOS 



The formal solution to Eq. ( 15 .2 1 1 is 

v(7) = v(0) exp 
2tc 2 



+ dl(2-z A ) 
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dl 



z Adl = — 



(5.22) 



With the notation e ren = e(/), e = e(0), we find 
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(5.23a) 



where g A (E ren ) — g A (e) is the solution to 
2k 2 



C) 2 A A (e) 



dgA^A- 



(5.23b) 



The renormalized energy scale E ren = e(/) is the energy scale 
at which we cannot neglect renormalization effects from the 
energy perturbation on the variances of the random potentials 
anymore. Beyond the energy scale e ren , the RG analysis en- 
coded by the Callan-Symanzik equation d5 .2 1 i breaks down. 
We take e ren as given, but we allow I to be tuned, or, equiva- 
lently, we take the bare energy e = e(0) to be tunable. How 
large the rescaling scale I is thus dictates how close the bare 
energy scale e is to the band center. The infrared limit / — > °° 
gives the DOS arbitrarily close to the band center. 

Equation (I5.23a> implies that the DOS is characterized by 
two regimes each of which can itself be divided into two sub- 
regimes. When 



8a(1)-8a(0) 
8a(0) 



«1, 



(5.24) 



the scale dependen ce of g A can be neglected in the Callan- 
Symanzik equation J5.21b> . We then recover the DOS for the 
MDH model J4.38i with za given by Eq. d5.21al . whereby the 
freezing transition can be observed as a function of the bare 
energy e depending on whether g A (e) i^ 2 > 2% or g A (e) i^ 2 < 271. 
When 



t (0-gA(0) 
*a(0) 



>1, 



(5.25) 



the scale dependence of g A in the Callan-Symanzik equation 
d5.21bl > is dominant. In the infrared limit / — > °°, i.e., arbitrar- 
ily close to the band center, the integral on the right-hand side 
of Eq. ( I5.23bl > is always dominated by the contribution 
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After solving Eq. d5 .23bl > for the leading dependence of g A on 
e ren , this gives the DOS 
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where c is some nonuniversal positive constant. Observation 
of the Gade singularity 



v(e) 
v(e ren ) 
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exp 
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(5.28) 



demands fine-tuning of bare and renormalized g A 's so as to 
avoid the freezing transition in Eq. I5.21at . At best, the Gade 
singular DOS can be observed as a crossover in a finite win- 
dow of energy. 



VI. CONCLUSIONS 

Two recent theoretical works have addressed the issue of 
the distribution of the DOS in the HWK model. Guruswamy 
et al. in Ref. |^ considered a narrow distribution of disorder 
at the microscopic level for which the use of field-theoretical 
methods at longer length scales is justified. They computed 
the DOS within a SUSY field theory and found that it dis- 
plays the Gade singularity £~'exp(— c| lne| x ' 2 ) in the neigh- 
borhood of the band center. Motrunich et al. in Ref. yCJ by- 
passing a purely field-theoretical approach, came to the con- 
clusion that the DOS is given by e _1 exp (—c\ lne| 2//3 ) instead 
of the Gade singularity. Our field-theoretical analysis agrees 
with the prediction of Motrunich et al. 

We find that the RG analysis performed on the DOS by 
Guruswamy et al. is not consistent as it overlooked the ex- 
istence of infinitely many operators with negative anomalous 
scaling dimensions. The physical interpretation of these op- 
erators is that the LDOS is a very broadly distributed random 
variable close to the band center. The situation here shares 
many similarities with the scaling with system size of the dis- 
tribution of the conductance in 2D mesoscopic samples which 
is known to broaden with increasing system size (i.e., by mov- 
ing away from the diffusive regime)..- 7 A parallel with quasi- 
one-dimensional wires is particularly instructive. In quasi- 
one-dimensions, it is known that the infinite set of scaling 
equations for the moments of the conductance can be traded 
for a Fokker-Planck equation obeyed by the underlying dis- 
tribution of the conductance^ The advantage of the Fokker- 
Planck equation over the scaling equations for the moments is 
that the Fokker-Planck equation can be solved in the localized 
regime (i.e., when the length of the sample diverges). It is then 
possible to show that the typical conductance is self-averaging 
(at least in the standard universality classes of Anderson local- 
ization) in the localized regime^2, Similarly, all moments of 
the LDOS are here strongly coupled through infinitely many 
RG equations encoding the dependence on energy arbitrarily 
close to the band center. The Gade singularity results from 
an inconsistent truncation of this infinite set of coupled equa- 
tions to a finite subset. As with the quasi-one-dimensional 
conductance, it is possible to trade the infinite set of scaling 
equations for the DOS in favor of a second-order partial dif- 
ferential equation obeyed by the underlying probability distri- 
bution for the normalization of zero modes from which one 
can extract the dependence on energy of the DOS. As with 



the quasi-one-dimensional typical conductance in the standard 
universality classes of Anderson localization, the DOS is well 
behaved statistically; i.e., it is self-averaging in the infrared 
limit. 

A classification of random Hamiltonians with short-range- 
correlated disorder in terms of intrinsic symmetries in con- 
junction with scaling concepts has emerged, since its incep- 
tion in the late 1970's, as perhaps one of the deepest in- 
sight into the physics of Anderson localization. The HWK 
model is invariant under time reversal and spin rotation. Those 
are the standard symmetries in Anderson localization. The 
HWK model also preserves a bipartite lattice symmetry, mak- 
ing it a member of the chiral orthogonal universality class. 
The HWK model is nevertheless special since the disorder is 
parametrized (in a way consistent with rescaling of the lat- 
tice spacing) by only two independent variances as opposed to 
four for a generic member of the chiral orthogonal universal- 
ity class. Since the works of Gade and Wegner, it is believed 
that the DOS in the chiral orthogonal universality class dis- 
plays the Gade singularity. If universality is indeed the only 
ingredient that determines the scaling exponent v and K for 
the singular density 



|er u exp(-c|lne| K ) 



(6.1) 



in the chiral orthogonal universality class, we must conclude 
that the RG analysis made by Gade in Ref. @ suffers from the 
same deficiencies as the RG analysis of Guruswamy et al. 

This possibility is perhaps not so surprising if we take the 
existence of infinitely many operators with negative anoma- 
lous scaling dimensions as the "smoking gun" for broadly 
distributed random variables. There have been indications 
since the early 1980's that replicated NLSM might support 
infinitely many operators with negative anomalous scaling 
dimensions^ A perhaps stronger clue comes from the in- 
terplay between multifractality of zero modes and the broad 
distribution of the LDOS in the HWK model. As we have 
shown explicitly in this paper, essentially the same spec- 
trum of negative scaling dimensions controls the multifractal 
spectrum of zero modes in the MDH model and the Callan- 
Symanzik equation for the DOS near the band center in the 
HWK model. Fal'ko and Efetov have searched for and found 
in the diffusive regime of 2D disordered metals the signature 
of multifractal states^i There are states, dubbed prelocalized 
states, that are characterized by very large amplitude fluctua- 
tions whose statistics are well described (in a finite window of 
length scales) by a log-normal distribution, i.e., by the multi- 
fractal spectrum of zero modes for a Dirac spinor subjected to 
a random vector potentiali22i&I It could very well be that the 
gl(l 1 1) current algebra at the heart of the freezing transition in 
the multifractal spectrum of zero modes in the MDH model 
and for the DOS of the HWK model, although broken by all 
the additional degrees of freedom encoded by NLSM in the 
standard or chiral universality classes, is sufficiently robust 
to reveal itself in those latter models in connection with pre- 
localized states and broadly distributed LDOS, respectively. 
To put it differently, the critical statistics of wave functions 
might still be log normal even though the wave functions have 
a much more complicated structure than a single exponential 
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of a field with logarithmic spatial correlations. To add credit 
to this speculation, we know of one example where critical 
modes have a much more complex structure than that of a 
single exponential of a two-dimensional scalar field (a Gaus- 
sian random surface in the terminology of Motrunich et al.) 
although their multifractal spectrum is still of the Gaussian 
random surface type: N flavors of Dirac fermions subjected to 
a random SU(N) gauge potentiali^LS An interesting direction 
that might be worth exploring is the construction of variational 
methods that would generalize to the NLSM the one devised 
by Horovitz and Le Doussal in Ref. |42| to capture the non- 
perturbative physics of freezing in the 2D random phase XY 
model. 

Another intriguing question is whether multifractality of 
critical states at a mobility edge or at the critical energy in 
the quantum Hall plateau transition necessarily reveals itself 
in a broad distribution of the LDOS. If the answer is positive 
as is suggested by numerical simulations in Refs.l63l and l64l 
the LDOS might be more interesting than commonly thought. 
Indeed, it has been proposed that the freezing transition in 
the multifractal spectrum of zero modes in the HWK model 
is generic to multifractal wave functions at a mobility edge 
and, in particular, also applies to critical states in the plateau 
transition^ The corollary would then be that a counterpart to 
this freezing transition characterizes the statistical distribution 
of the LDOS — say, that typical and average LDOS disagree. 
So far, relatively little attention has been paid to the distribu- 
tion of the LDOS in the plateau transition since the calculation 
by Wegner of the DOS (Ref. l6o) in spite of great theoreti- 
cal efforts to identi fy th e critical theory governing the plateau 
transition (see Refs. l67l andl68l for the latest endeavors). 

Note added in proof. Within the SUSY representation of the 
Gade model (a NLSM) introduced by Guruswamy et al. it can 
be seen that it is the so-called Gade term which is, ultimately, 
responsible for the freezing transition of the DOS. The Gade 
term in the Lagrangian of the NLSM results from the fact that 
the superdeterminant of the supermatrix T entering the NLSM 
is not constrained to some fixed numerical value as is the case 
in all other universality classes. The superdeterminant of T 
is thus a dynamical variable of its own whose fluctuations are 
controlled by the Gade term (3 ;U lnSdetr) 2 . A similar role is 
played by the Gade term in the bosonic and fermionic repli- 
cated NLSM representing the chiral universality classes. 
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APPENDIX A: THE PROCESS OF COMPLETE 
ANNIHILATION 

In this appendix, we give the full OPE's obeyed by the com- 
posite operators that control the DOS in the MDH model of 
Sec. lIVI Although this paper deals with two flavors of Dirac 
fermions subjected to an imaginary random vector potential 
and a complex-valued random mass, we also consider the 
case of one flavor subjected to a real-valued random vector 
potential^ Indeed, as we have argued, the problem of the 
DOS in the two-flavor problem reduces to a large extent to 
the one-flavor problem. Moreover, the one-flavor problem is 
interesting in its own right as it is related to the plateau transi- 
tion in the integer quantum Hall effect? and to the 2D random 
phase XY modeli 27 i 35 i 37 i 41 



1. Case of one flavor 

The problem of a single flavor of Dirac fermions coupled to 
a real-valued random vector potential can be shown to reduce 
to the action^ 

2 8 A J 

+ j -^«5y_+v$-3v+) 

+ / — (p*_8p_ + p;ap + ). (ad 

Here, 4>i is the transversal component of the random vector 
potential, with variance g A . For simplicity, the longitudinal 
part <t>2 and the Jacobian that arise when trading the vector 
potential by its transversal and longitudinal components (con- 
formal sector with net central charge —1) have been dropped 
from the action as we shall only seek the scaling dimensions 
of gauge invariant composite operators. Action \A\\ is critical 
for arbitrary value of g A as shown in Ref. Band defines a CFT 
with gl(l 1 1) symmetry as shown in Ref.l37l 

It was shown in Ref. |^| that there exists an infinite set of 
operators with negative anomalous scaling dimensions. These 
operators are composite operator of the form 

O m ,n = {<F + ) m+ {* + ) n+ {?-) m - («-)"- 

x exp [2 [m+ +«+ — m_ — «_)<t>i] , 

(A2) 

where f± and r B± are defined as 

S+:=P*_P + , S_:=p;p_. (A3) 

Their scaling dimensions are given by 

■W = ( m + - m - ) 2 + \n+ - «- 1 

-^[(m + -m_) + (n+-n_)] 2 , (A4) 

K 
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respectively. We thus see that for any value of the disorder 
strength g A , there are infinitely many operators with negative 
scaling dimensions as (n+ — n_) 2 dominates over |n + — «_| 
for large \n+ —n~\. Of all the composite operators {o m ,n}, 
we focus on the subset spanned by 

(<B ± ) N e ±2N *2, (f ± )(v ± ) N - l e ±7 - N *2 7 (A5) 

where N isa positive integer. For any given N, these two oper- 
ators are the most relevant and share the same scaling dimen- 
sion, forming a two-dimensional representation of gl(l 1 1). We 
then construct gl(l 1 1) sing lets21 

A_:Ar:=_ 3vr (®- +Nyr - rB - )> < A6 > 

whose scaling dimensions read 

x±- N = x N =N— —N 2 . (A7) 
Jt 



Operators (IA6> are always induced in a RG analysis of the 
DOS if one is after localization properties or in a RG analysis 
of the relevance of charge 2 vortices if one is after the fate of 
the Kosterlitz-Thouless transition in the 2D random phase XY 
model. Indeed, they are induced through repeated OPE's of 
the operator 



D=e +2ct,1 (2i + + J + )+ e - 2 * 1 (2!_ + J_) (A8) 



with itself. Hence, any RG analysis of the DOS or of the 
stability of the spin-wave phase in the 2D random phase XY 
model is controlled by the OPE's within the set {.As-jv}- These 
OPE's can be classified into three different processes which 
we call "fusion," "partial annihilation," and "complete annihi- 
lation" and are given by 
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^ + (z,z)^_(w,w) = (-) N -\z-w\ XN+ - N - ~ Xn +~ Xn - ^ N+ N ~ 1 ^jx + ; N+ - N _(w,w) + ---, 

0<N_<N + eN, (A9b) 

A + ;N + {. Z ->Z) A -\N-{W,W) = {-f \z~w\ 2 ~ 2x "{ + {N -l)0 A (w,w)+NO M (w,w) 

-2Nd<t>i(w) [-f(w)+J(w)] + 2Nd<i>i (w) [-J'(w)+J(w)] } + •••, 

0<N_=N + =NeN, (A9c) 



respectively. These OPE's are the central result of this section. 
We have used the shorthand notation 

o A : = -« ¥ _+p*-P-)(v;¥ + +p;p + ) 

= -(/'-/) (/'-/) (A10) 
for the random vector potential perturbation and 

o M := -( ¥ ; ¥ _+p;p_)« ¥+ +p*-P+) 

= + (//- // -G-G+ + G+G ) 

(All) 

for the random mass perturbation. The gl(l|l) currents are 
defined by 

/:=\|/_\|/*_, /:=p_p*, 
G+:=p*_ ¥ _, G-:=$_\f_, (A12) 



with similar definitions in the antiholomorphic sectors for 
/,/',G+,G_. Operators Oa and Om are induced by disor- 
der averaging over a Gaussian-distributed real-valued random 
vector potential and random mass, respectively. We see that 
the process of complete annihilation induces a renormaliza- 
tion of gA through the presence of the terms 3<t>i (—/' + /), 
34>i (— / +f), and A on the right-hand side of Eq. jA9cl We 
can close the OPE's ( IA9> by inclusion of the OPE's {s = ±) 

/ -x / N(N-l) . , 

Om{Z,Z)As;N(W,W) = 7^-M. s - t N(W,W) H , 

\Z — W\ A 

(A13a) 

-2 

M (z,z)0 M (w,w) = [T Oa(w,w)H . 

\Z — W\ z 

(A13b) 



25 



It is consistent to neglect the process of complete annihilation 
jA9^ if one neglects OPE's jA13> . This is the approximation 
that we make to calculate the DOS arbitrarily close to the band 
center. 



2. Case of two flavors 

We now return to the OPE's of the MDH model in Eq. 
(14. 1 3i which we supplement by specifying the process of 



Here, 

oa = - I «_i|/ a _+p:_p fl _)(^ + i|/ 6+ +p; + p 6+ ) 

a,b=l 
2 

52 [ ~^aa^bb + ^ aJbb + ^aa^bb ^ a<Jbb\ 
a,b=l 

(A15) 

is the random vector potential perturbation for the two-flavor 
case, and Oo is an operator generated by taking the OPE 
J4. + -\ x whose explicit form is given by 

Oo = 

+•^11^22 + ^12^21 + ^21^12 +^22-^ll 

— Gi 1+G22+ — G\2+Gi\+ + G21+G12+ + G22+G11+ 

— G1I-G22- — G12-G2I- + G2I-G12- + G22-Gn_ 

—Jnfzz+Jn-hi +H\Jvi — hih\- (A 16) 

Observe that the multiplicative factor (— ) N present in Eq. 
dA9> does not appear here. This is so because we find op- 
erators with negative scaling dimensions in the gauge-variant 
sector, i.e., through "Cooper-like" terms, when the random 
vector potential is imaginary. As was the case for a real- 
valued random vector potential, a renormalization of gA is 
encoded by the complete fusion process through the contribu- 
tions 34>2l„ (-J'aa+Jaa) , 3*2 L, (-J' aa + J aa ), and A . For 



complete annihilation 



(A14a) 



(A14b) 



(A14c) 

I 

the MDH model, we neglect this effect in the same way as we 
neglect the renormalization of g A induced by the source term 
needed to compute the DOS. This approximation should be 
good if one is after the DOS asymptotically close to the band 
center but fails if one is interested in the DOS at large energy 
scales. For the HWK model, we neglect complete annihilation 
relative to the renormalization of g A induced by g M . 

APPENDIX B: BOSONIZATION OF THE GHOST SYSTEM 

1. Bosonization 

The purpose of this section is to review the steps of Ref.l27l 
by which the SUSY action 

f d 2 r 

So = J — (vi3y_ +Y+3y+) 

+ |^(p*_ap_ + p;ap + ), (Bi) 

where are fermionic spinors and p*,p are bosonic 

spinors, can be "bosonized," i.e., represented by a pair of 
independent scalar fields (j) and (j)' and a pair of independent 
fermionic fields %* and %. The fermionic sector in the SUSY 
action JBli describes a CFT with central charge c = +1. The 
bosonic sector in the SUSY action JB U describes a CFT with 
central charge c = — 1 . Taken together, both sectors describe 



( N + N' \ 

\z - w\ N+N ' N n ' ( N J *«ft+N> (w, w) + ■ ■ ■ , 

0<tf,iV'eN, s = ±, 

\Z - w\ X » + ~»- - XN + - XN - ( N+ N ~ 1 ) * + ; N+ - N _ (w, w) + ■ ■ ■ , 

0<N- <N+en, 
|z-w| 2 ~ 2v ~{ - (N— \)o a (w,w)+NOq{w,w) 

2 2 
+2Nd® 2 {w) £ [-?aa(*>)+Jaa(w)\ +2Nd® 2 (w) £ [-j' aa (w) + J aa ( W )] } 



a=l a=l 

0<N+=N-=Ne'N. 
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a CFT with vanishing central charge. The OPE for the spinors 
are 

¥ _(z) ¥ *_(0) - <(z)V-(0) = - + •••, 



V+(z)v;(0) = ¥* + (z) ¥+ (0) = r 

z 



P_(z)P*_(0) = -P* (z)P_(0) = - 

Z 

p + (f)p;(o) = -p;(f)p + (o) = i- 



After defining the dimension- 1 currents 



(B2) 



/:= ¥ _ ¥ *_, /:=P_P*, 

/:= ¥+¥ ;, /:=P+P;, 

G + :=p*_ ¥ _, G_:=p_ ¥ *_, 

G + :=p; ¥+ , G_:=p +¥ ;, (B3) 

we see that they generate the gl(l 1 1) super current algebra 

J(z)J(w) 
J'(z)f(w) 



(z~w) 2 

-1 

(z-w) 2 

1 



J(z)G±{w) = ± G±(w) 

z — w 

j'{z)G±{w) = ±— ^G±(w)- 
z — w 

G T (z)G ± (w) = 

-1 



+ [7(w)-7'(w)] + . 

z — w 



(B4) 



The idea behind bosonization is to represent the SUSY ac- 
tion ( IB 1 b in terms of a new set of fields so as to preserve the 
OPE's J B 2b and ( IB4I . For example, it has been known since 
the early 1970's that the currents for the fermionic c = 1 sector 
can be mimicked by a scalar field (j)(z,z) through the identifi- 
cations 



J /=-(i5$), 



(B5) 



which preserve the J — J OPE in Eq. JB4> provided (j)(z,z) 
obeys the OPE 



<j)(z,z)(j)(0) = -lnzz+. 



(B6) 



In turn, fermionic spinors can be expressed by the Mandel- 
stam formulas 



>|fl = (M*) +1 a- 1 / 2 e + * l4 - > 



(B7) 



whereby (j)(z,z) is decomposed into a holomorphic and anti- 
holomorphic part according to (j)(z,z) = (p_ (z) +cp + (z). Here, 
an arbitrary complex number u has been introduced, as it is 
not fixed by the OPE's dB2l . 

Although the p* — p subsystem is already bosonic, it can 
be advantageous to trade these spinors for another scalar field 
<j)'(z,z) = 9-(z) +cp' + (z)i^ As in the fermionic case, identifi- 
cation 



(B8) 



is compatible with the /' — /' OPE in Eq. (IB4> provided (j)' 
satisfies 



<j)'(z,z)(j)'(0) = +lnzz- 



(B9) 



Observe that this OPE differs from the one in the (j) sector by a 
sign. Since the p* — p sector can be thought of as a CFT with 
central charge c = — 1 whereas it is natural to assign the central 
charge c = +1 to the (j)' sector, a conformal sector with central 
charge c = —2 must be supplemented to the 0' sector. The 
need for an additional conformal sector can also be seen when 
tentatively writing down Mandelstam formulas for the bosonic 
spinors P±, P± of the form iBli with the substitution of cp ± by 
(p' ± as the exponentials exp(±/(p' ± ) obey Fermi statistics. The 
correct bosonic statistics in the P* — P sector is implemented 
after the introduction of two pairs of independent fermions 
ri± and £± in the holomorphic and antiholomorphic sectors, 
respectively, through the bosonization rules 



P_ = v +1 
p*_ = y-' 



P + ^(v*)+ 1 a +1 /2 e -V +11+ , 

p; = (v*)- 1 a +1 / 2 e +i <^ + . 



(BIO) 



As before, the arbitrary complex number v has been intro- 
duced, as it is not fixed by the OPE's iB2\ . whereas T|_, 
T| + , and t,+ are anticommuting fields with the conformal 
weights (1,0), (0,0), (0,1), and (0,0), respectively. Their 
OPE's are thus given by 



S-(Z)T1_(0) = Tl_( Z )$_(0) = i 

^(z)T! + (0) = iu(z)UO) = 4 



(Bll) 



Finally, the action dBll has the bosonized representation 
d 2 r 



s = I ^(ew-W) 

d 2 r 



(B12) 



Bosonization of the p* — p sector is not unique — Instead of 
the r| — i; sector we could have chosen to introduce (symplec- 
tic) fermions through the action 



d 2 r 



(B13) 
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the OPE's 

X*(z,z)%(0) = -%(z,z)X*(0) = -buz- 
and the identifications 



n-(z) = '3x*(z). 
n+(z) = ih*{z), 



^-(z) 



'3z(z), 
'3x(z)- 



(B14) 



(B15) 



[The minus sign in Eq. JB 1 3i is necessary to get the minus 
sign on the right-hand side of the OPE IB 141 with the con- 
vention we are using for Grassmann integration.] With the 
Mandelstam formulas, fermionic currents are bosonized as 



G 
G 

G+ 
G+ 



(uv) \e +i ^'--^-\idx*) 

(uv)* le-W+'^iidx*) 

(uv)' 1 le-^'-^iidx ) 

[(uv)- 1 ]* \e +i ^-^(id% ) 



(B16) 



This bosonization rule, together with Eqs. ( IB5i and ( IB8I . is 
compatible with the gl(l 1 1) super current algebra OPE's (IB4> 
for any value of u, v £ C. 

The ambiguity in the choice of the phase of (uv) when 
bosonizing fermionic currents is, however, lifted in the 
presence of gl(l|l) current-current interactions by imposing 
boundary conditions on the bosonized fields that maintain 
SUSY. To see this, consider the action [recall Eq. (!5.5bH 

c c 8A fd 2 r g M fd 2 r 

S = S0 -^J— OA{r) -2^Jl^ OM{r) > 

(B17) 

which describes action JB 12i perturbed by the random vector 
potential Oa and the random mass Om defined in Eqs. ( IA10I 
and JA1 1> . respectively. Although S is not conformally in- 
variant it is still GL(1|1) symmetric. Consequently, the cor- 
relator (Om(z,z)Om(w,w)}, say, vanishes for any gA and gM- 
Calculation of this correlation function can be done nonper- 
turbatively as a function of gA and gM along the lines of Ref. 
I27L The condition that this correlator vanishes uniquely speci- 
fies the boundary condition obeyed by the scalar fields and, in 
turn, the phase of (uv). In this paper we choose the boundary 
condition for which the bosonized fields vanish at infinity in 
the CFT defined by 



d 2 r 

— (dx*dx+dx*dx). 



(B18a) 



delstam formulas become (with \uv\ = 1, u = v, say) 

¥ _ = a -l/2 e -m/4+'-<P- 7 

^|f=a- 1/2 e + ' ,l < /4 -'' ( f^ 
i|/ + = cr 1/2 e+" t / 4 - , ' < f > +, 
\ V * + = a- l/2 e-^ 4+i( f + , 



(B18b) 



and 



P- 



+ 1/2 +ra/4+,V_ 



^_ = a +l / 2 e-' K / 4 -^'-(idx), 
^, + = a +x l 2 e- in l A - i ^(idx*) 1 
p; = a +1 / 2 e + " t / 4+, ' < P'+(/3x). 



(B18c) 



2. Scaling dimensions of composite operators made of bosonic 
spinors 

The x sector plays an essential role as it ensures that the 
P sector obeys a bosonic algebra. It also guarantees that the 
central charge of action (IB18a> vanishes and gives compos- 
ite operators from the (3 sector their correct scaling dimen- 
sions. To illustrate this last point we calculate the scaling di- 
mension of the composite operator (P+P_ )". Taking repeated 
Wick contractions, one verifies that the scaling dimension of 
(P^_P_)" is n. On the other hand, bosonization rules JB18I > 

give (P^P_)" = a +n i n e +ini >' (%*)" (idx)". Since e +in6 >' con- 
tributes the negative scaling dimension — n 2 in view of the 
OPE obeyed in the (])' sector, (idx*) (f^X) must contribute 
the anomalous scaling dimension n 2 + n. That this is indeed 
so follows from noting that is generated by repeated 

OPE of (idx*) with itself , 



(idX*T(z) 



1 



1 1 



:3* - W) 



(n-1)! 1!0! 
xd"- 2 (idx*)---d l (idx*)d°(idx*):(z). 

(B19) 



[The notation :(■••): refers to normal ordering.] Since the 
conformal weight of (idx*) i s (1,0), the conformal weight of 
(idx*)" reads 



(n+(«-l)+...+ l,0)=('^±^,0 



This amounts to the choice uv/\uv\ = —i for which the Man- 



(B20) 



Accounting for the contribution from the antiholomorphic part 
(idx) yields the anomalous scaling dimension n(n + 1) for 

(idx*)" fix)" ■ 
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